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Preface 



The twentieth century witnessed the development of both gravitation and cosmology as 
modern scientific disciplines subjected to observations. These observations have been per- 
formed through terrestrial particle detection devices, telescopes and sateUites that allow 
to verify theoretical predictions and to rule out proposed theoretical models. With the 
turning of the new century, called to be the century of precision cosmology, new pers- 
pectives have been unveiled with recent experiments such as WMAP, PLANCK or SDSS. 
These last experiments are able to determine with higher and higher accuracy the features 
of the Cosmic Microwave Background (CMB), the distribution of large scale structures 
and the fundamental cosmological parameters which describe our universe on the largest 
scales. Despite the improvements in the observational side, a fundamental gravitational 
theory, which is renormalizable from a quantum field theory point of view and applicable 
to arbitrary scales, from micro-gravitational tests, passing through solar system tests, to 
cosmological scales, is still lacking. 

General relativity, in spite of being the most successful gravitational theory in the last 
one hundred years, has left some of these problems without satisfactory answer. Although 
within the string theory paradigm it would be possible to find a consistent quantum theory 
of gravity, this is not the case of general relativity which turns out to be nonrenormalizablc 
as a perturbative field theory. Moreover, if this theory is used to construct the standard 
cosmological model, where the fiuid content is given by standard matter and radiation, 
it cannot account for the observed accelerated expansion of the universe on sufficiently 
large scales. In fact, it needs to be supplemented by some dark energy contribution to ac- 
commodate this accelerated regime. On the other hand, general relativity with gravitating 
luminous matter cannot account either for the observed rotation curves of galaxies. A dark 
matter contribution needs to be introduced to reconcile data with theoretical predictions 
within this paradigm. 

Instead of adding new elements in the cosmological content, which try to accommo- 
date observations with general relativity, those problems might show that the theoretical 
framework in cosmology should be enlarged by alternative gravity theories. This thesis 
will try to contribute to the understanding of those still open issues by considering two 
recently proposed alternative and complementary theories to general relativity. We shall 
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consider some relevant aspects of those models related to recent experimental results. 

The present work is organized in the way that follows: First, we will briefly introduce 
in Chapter 1 some modified gravity theories and their corresponding formalisms. In this 
chapter special attention will be paid to f{R) gravities by summarizing the main features 
of this paradigm in the metric formalism. Then some geometrical results for f{R) theories 
and both cosmological and gravitational constraints usually imposed over such functions 
will be provided. Other alternative modified gravities, the brane worlds, will then be 
reviewed. Here we shall introduce both the notion of brane excitations, the branons, and 
some topologically nontrivial solutions, the branc-skyrmions. We shall finish the chapter 
by providing some insight about the possibility of mini black holes detection in the Large 
Hadron Collider (LHC) as a signature for the validity of these modified gravity theories. 

The second chapter will deal with f{R) theories which try to provide a cosmological 
acceleration mechanism with no need for introducing any extra dark energy contribution in 
the cosmological components. To do so, we shall use some reconstruction procedures which 
start either from a given solution of the cosmological scale factor for an homogeneous and 
isotropic metric or from an effective equation of state. In particular, those /(it!) functions 
able to mimic Einstein-Hilbert plus cosmological constant solutions will be obtained. In 
this realm, f{R) theories will be shown to be able to mimic the cosmological evolution 
generated by any perfect fluid with constant equation of state. 

Then the third chapter will be devoted to the computation of cosmological perturba- 
tions for f{R) theories. Since in Chapter 2 the modified Einstein equations will have been 
studied as background equations, it is quite natural when modifying general relativity by 
f{R) models, to ask about the first order perturbed equations for these theories and what 
consequences in the growing of these perturbations may appear. This is the leitmotiv in 
this chapter. Throughout it, special attention will be paid to the possibility of obtaining 
a completely general differential equation for the evolution of perturbations and its par- 
ticularization for the so-called sub-Hubble scales will be explicitly shown. The mentioned 
differential equation in those scales will be shown to be very useful to understand the 
regime validity of some approximations widely accepted in the literature and to rule out 
that some proposed f{R) models could be cosmologically viable. 

The introduction of modified gravity theories, with or without extra dimensions, may 
lead to the existence of new solutions with respect to those of general relativity. In that 
sense, the research about spherically symmetric solutions is of particular interest. For 
instance, it may shed some light on the number of extra dimensions, the fundamental scale 
of gravity or the required restrictions to be imposed over the parameters of those theories. 
The possible detection of mini black holes at the LHC in the coming years will be a turning 
point to discover certain properties of the underlying gravity theory. For this reason, 
chapters 4 and 5 will be devoted to the study of spherical solutions in extra dimensions 
theories. In particular, spherically symmetric and static black-hole solutions coming from 
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f{R) theories in an arbitrary number of dimensions will be studied in Chapter 4 whereas 
Chapter 5 will be focused on studying a particular topologically nontrivial solution with 
spherical symmetry appearing in brane- world models - different from the well-known black- 
hole solutions - the so-called brane-skyrmions. 

Hence, in the fourth chapter we shall focus on the study of black holes in f{R) gravity 
theories in an arbitrary number of dimensions. Wc shall concentrate on the existence of 
black-hole solutions and we shall study which will be their inherited or different features 
with respect to those in general relativity. With this purpose we shall study constant cur- 
vature solutions for f{R) theories as well as perturbative solutions around the standard 
Schwarzschild-anti-dc Sitter geometry. An important part of this chapter will be then de- 
voted to the thermodynamics of Schwarzschild-anti-de Sitter black holes in /(/?) theories. 
This research will prove that for f{R) gravities, there exists a thermodynamical viability 
condition which is related to one of the conditions which ensure gravitational viability for 
f{R) models. 

In the fifth chapter we will thoroughly study other kind of spherically symmetric 
solutions in brane- world theories that are not black holes. These solutions, the brane- 
skyrmions, are topologically nontrivial configurations arising in the presence of these extra 
dimensions theories. In this context, the recent claim of detection of an unexpected feature 
in the CMB, referred to as the cold spot, will be explained as a topological defect on the 
brane. After performing some calculations, it will be shown that results obtained are in 
complete agreement with those in the literature that tried to explain that cold spot as 
a texture of a non-linear sigma model. The physical interpretation of these results and 
future prospects will finish this chapter. 

At the end of each chapter, we shall include the corresponding conclusions. These 

conclusions arc summarized all together in the sixth chapter, which is followed by an 
appendix where more detailed formulae for the calculations performed in the third chapter 
are shown. 
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Chapter 1 

Introduction to modified gravity 
theories 

1.1 Motivation 

From its very beginning, it was questioned whether general relativity (GR) was the unique 
correct theory among other theories for gravitation. Thus for instance Weyl in [1] and Ed- 
dington in [2] included higher order invariants in the gravitational action. Those attempts 
were neither experimentally nor theoretically motivated, but it was soon proved that the 
Einstein-Hilbert (EH) action was not renormalizable and therefore could not be conven- 
tionally quantized. In fact, this action needs to be supplemented by higher order terms in 
order for the resultant theory to be one-loop renormalizable [3, 4]. More recent research 
has shown that when quantum loop corrections in field theory or higher order corrections 
in the low energy string dynamics are considered, the effective low energy gravitational 
action includes higher order curvature invariants [5, 6, 7]. 

Such results encouraged the interest in higher order gravity theories, i.e., modifications 
of the EH gravitational action which include higher order curvature invariants. Nonethe- 
less, those new added contributions were thought to be relevant only in very strong gravity 
regimes, such as at scales close to the Planck scale and therefore in the early universe or 
near black hole singularities. However, these corrections were not expected to affect gravi- 
tational phenomenology neither at low curvature nor at low energy regimes, and therefore 
they were assumed to be negligible at large scales such as those involved in the late universe 
evolution. 

Very recent evidence coming from both astrophysics and cosmology have revealed the 
unexpected accelerated expansion of the universe. Different data from type la supernovae 
(SNIa) surveys [8, 9, 10], large structure formation and delicate measurements of the CMB 
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anisotropies, particularly those from the Wilkinson Microwave Anisotropy Probe (WMAP) 
[11], have concluded that our universe is expanding at an increasing rate. This fact sets 
the very urgent problem of finding the cause for this speed-up since standard GR with 
ordinary matter and radiation is not able to do so. Usual explanations for this fact have 
been categorized to belong to one of the following three classes: 

1. The first type of explanations reconciles this acceleration with GR by invoking a 

strange cosmic fiuid, dark energy (DE) (see [12] and references therein), with a state 
equation relating its pressure and energy density in the following way 

-Pde=^^dePde (1-1) 

where cjde < —1/3 is required to provide acceleration in the usual Einstein equa- 
tions as is described in Section 2.2. This state equation shows that the DE fiuid 
has a large negative pressure. For the particular case cude = — 1 , this fiuid behaves 
just as a cosmological constant A . Within this approach of DE in the form of a 
cosmological constant, recent data obtained by WMAP [11] provide the following 
cosmological content distribution: 4.6% corresponds to ordinary baryonic matter, 
22.7% to cold dark matter and 72.9% to DE. This is the so-called concordance or 
A -Gold Dark Matter model (AGDM) which is supplemented with some infiation 
mechanism usually through some scalar field, the infiaton. The main problem of 
this kind of description is that the fitted A value seems to be about 55 orders of 
magnitude smaller than the expected vacuum energy of matter fields, this is the 
so-called cosmological constant problem. Prom a more philosophical point of view, 
the DE description also presents the so-called coincidence problem. This problem 
wonders why the DE and matter densities are so close in order of magnitudes pre- 
cisely in these days, i.e. in the present cosmological era, even though for both the 
cosmological past and future that is not the case. This kind of problems comes to 
claim that the AGDM model could be regarded as an empirical fit to data with a 
poorly motivated gravitational theory behind and therefore, it should be considered 
as a phenomenological approach of the underlying correct cosmological theory. 

2. The second type of explanations consider a dynamical DE by introducing a new 
scalar field. They are the so-called quintessence theories. The theories which intro- 
duce an extra scalar field in the gravitational sector of the action are usually referred 
to as scalar-tensor theories. Some very interesting subcases of such theories are the 
so-called Brans-Dicke theories which are going to be explained in detail in the Section 
1.4. 

3. Finally the third one consists of trying to explain the cosmic acceleration as a 
consequence of new gravitational physics [13, 14]. For instance, modifications 
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to the EH gravitational action have been widely considered in the hterature 
[15, 16, 17, 18, 19, 20]. More recently, vector-tensor theories of gravity and the 
electromagnetic field itself have also been proposed as compelling DE candidates 
[21]. 

Some of those theories add higher or lower powers of the scalar curvature, the Rie- 
mann and Ricci tensors or their derivatives [22]. Lovelock theories and f{R) gravity 
theories are some examples of these attempts. In recent years, some f{R) proposals 
have even tried to reconcile dark matter through a gravitational sector modification 
[23] or to explain both the current cosmic speed-up and early inflation simultaneously 
[24]. The core of Chapter 2 will thoroughly deal with some attempts of f{R) theories 
to circumvent the necessity of introducing DE to explain the cosmic acceleration. 

On the other hand, other open issues in the Standard Model (SM) of elementary parti- 
cles, namely the hierarchy problem, could also be related to the fundamental gravity theory. 
Thus, this problem appears in the renormalization procedure in theories containing scalar 
fields. In such theories the renormalizcd scalar masses are expected to be given by the 
cut-off of the theory, i.e., the Planck scale. Therefore an extreme fine tuning is required in 
order to get the expected mass for scalars, in particular the Higgs mass. If on the contrary 
the fundamental scale of gravitation is close to the electroweak scale, the corresponding 
cut-off would be of the same order as the expected Higgs mass and an extreme fine tuning 
would not be required. 

With the aim of solving this problem, large extra dimensions theories have recently 
been considered. Unlike ancient Kaluza-Klein theories, with compactified Planck scale size 
extra dimensions, recent brane-world models may contain much larger extra dimensions. 
In order to avoid the presence of Kaluza-Klcin towers of copies of SM particles with similar 
masses, these models restrict SM particles to propagate on the brane, whereas only gravity 
can propagate in the whole bulk space. In this way, the fundamental gravity scale can be 
reduced to the electroweak scale and the gauge hierarchy problem is avoided. 

Brane-world (BW) theories may also explain the observed accelerated expansion of 
the universe [19] and as will be shown in Section 1.8, they present excitations which can 
produce weakly interacting massive particles (WIMPs), which are natural candidates for 
the observed dark matter [25]. Let us finally remark that such modified extra dimensions 
gravity theories, as will be explained in Section 1.10, may give rise to the production of 
black holes (BHs) of little size at the LHC whose eventual detection may give valuable 
information about the dimensionality of space-time. 

In the following sections of this chapter we shall deal with different aspects of the 
already mentioned modified gravity theories, both f{R) theories and brane-world theo- 
ries: in Section 1.2 we shall present some generalities about the formalism which will be 
used throughout the thesis. Thus in Section 1.3, the modified Einstein equations derived 
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from /(it!) theories in the metric formahsm will be presented. Then in Section 1.4, the 
equivalence of such theories with Brans-Dicke theories will be briefly sketched. Some geo- 
metrical results for f{R) gravities which were originally published in [26] will be presented 
in Section 1.5. They deal with constant curvature solutions and analytical conditions to 
reproduce Einstein's equations for the EH action with or without cosmological constant. 
Concerning BW theories, the main concepts for those models are analyzed in Section 1.7, 
while a study of BW excitations, called branons, is presented in Section 1.8. Topologi- 
cally nontrivial configurations, called skyrmions are introduced in Section 1.9 and some 
gravitational consequences of those theories will be summarized in the final Section 1.10. 



1.2 Generalities 

The gravitational action for GR in an arbitrary number of dimensions D is given by the 
so-called EH action 



Here, k = SttGd where Gd = M"^^ holds for the D -dimensional gravitational constant, 
with Md the gravitational fundamental scale, g is the metric determinant and R is the 
Ricci scalar defined from the metric tensor. 

With the aim of modifying the EH action, gravitational action for f{R) theories, con- 
sidered as generahzations of GR, may be written as 



Prom either actions given in equations (1.2) or (1.3), the field equations, giving rise to 
the so-called standard and modified Einstein equations respectively, can be derived by 
using different variational principles. Two such variational principles have been mainly 
considered in the literature: on the one hand, the standard metric formalism considers that 
the connection is metric dependent and therefore the only present fields in the gravitational 
sector are those coming from the metric tensor. On the other hand, there exists the 
so-called Palatini variational principle where metric and connection are assumed to be 
independent fields. In this case the action is varied with respect to both of them. Whereas 
for an action linear in R such as that in expression (1.2) both formalisms lead to the 
same field equations, this is no longer true for nonlinear gravity theories (see [27] for an 
exhaustive review on nonmetric formalisms). In this thesis, we shall restrict ourselves to 
the metric formalism. Por that purpose, we shall assume that the connection is the usual 
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r% ^ i^"^ {^ + ^-^) (1-4) 



Levi-Civita connection given by 

where, as in the rest of the work, Einstein's convention for imphcit summation is assumed. 

At this stage, let us point out that the convention to be used for the metric signature 
will be (+, — , — ) , i.e., positive sign for temporal coordinate whereas negative sign for 
spatial ones. With respect to the Riemann tensor definition, our conventions will be 

pM — I -per _ -pfi -pa (-t tC\ 

From expression (1.5), the corresponding Ricci tensor and scalar curvature are obtained 
straightforwardly and they read respectively as follows 

= ; R = R°a- (1.6) 

In addition to the already explained gravitational sector, the energy content may be 
introduced in the cosmological content through energy-momentum tensors, which will des- 
cribe the different components such as dust matter, radiation, dark matter, etc. which 
are present in the cosmological content of the universe. For each different type of fluid 
content ( a ), assumed from now on to behave as a perfect fluid, the corresponding energy- 
momentum tensor is given by 

T^S^ = (Pa + PaX^ u<;^^ - g^, (1.7) 

where Pa , Pa and are the pressure, energy density and 4- velocity of the a com- 

ponent respectively. Therefore the total energy-momentum tensor will be nothing but 

a 

for all possible fluid contributions. The most usual approach is to consider barotropic 
fluids where Pa = Pa{Pa) and very often the relation between these two quantities is 
linear through an equation of state 

Pa = UJaPa (1-9) 

where for instance Wq, = — 1,0,1/3 if cosmological constant, dust matter or radiation are 
the considered fluids respectively. DE fluids with constant equation of state are given by 
the condition cude < — 1/3 whereas phantom candidates for DE obey cjde < — 1 ■ In our 
approach to modify GR, to be rigorously implemented in Chapter 2, DE will appear as a 
modification of the gravitational sector itself so no DE component will be explicitly included 
in the content expressed by the summation (1.8). In this case the cosmic acceleration will 
be a consequence of the modification of the gravitational action by the presence of a f{R) 
term. Let us finish this section by mentioning that each fiuid component is assumed to be 
conserved separately since no interaction among fiuids is considered. This fact also implies 
the conservation of the total energy-momentum tensor straightforwardly. 
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1.3 Modified Einstein equations 

Now that the previous generahties have been presented, the modified Einstein equations in 
the metric formahsm for f{R) gravity theories may be found by performing variations of 
the gravitational action (1.3) with respect to the metric and equahng the result to minus 
the energy-momentum tensor times k providing the following equations: 

(1 + fR)R^. - liR + fiRM. + t>^Jr = -i^T^. (1-10) 

where Jr = df{R)/dR and 

V^, = V^V, - g^,n (1.11) 

with □ = Vq,V" and V is the usual covariant derivative. 
Taking the trace of the equation (1.10) we get: 

R{l + fR)-^{R+f{R)) + {l-D)DfR = -kT (1.12) 

which provides a differential relation between R and T unlike GR where this relation 
is just algebraic. An interesting point to stress at this stage is that in general, vacuum 
solutions, i.e. T^^, = , do not imply straightforwardly R — solutions. 

By computing the covariant derivative of (1.10), it is found that the l.h.s. of those 
equations vanishes identically, so the covariant derivative for the r.h.s. of equations (1.10) 
must obey the conservation equations 

v^n = (1.13) 

where this identity does not depend explicitly on f{R) but only on the energy-momentum 
tensor components and metric tensor elements. 

Two particular simple choices for f{R) may be considered in the equations (1.10): 

1. f{R) = , which allows to recover the standard Einstein equations without cosmo- 
logical constant, i.e., 

Gfj,u = R,_a, — -Rqi^u — —kT^i, (1-14) 
where the conservation equations (1.13) still hold. 



1.3. MODIFIED EINSTEIN EQUATIONS 11 



2. A second simple choice would be f{R) = — (-D — 2) Ajj . This choice allows to recover 
the standard Einstein equations in D dimensions with nonvanishing cosmological 
constant , i.e., 

1 D - 2 

R^u - ^Rg^u H — ^ — ^Dg^u = -i^T^u (1-15) 

where the particular choice of the Ad normalization will be explained below. Let 
us note that the equations (1.13) again hold. Notice that in this case the new piece 
in the previous equation (1.15) proportional to A^ can be moved to the r.h.s. and 
then an energy-momentum tensor (Ta^,)^^ can be defined as follows 

In this case, both density and pressure from the cosmological constant contribution 
may be written for any number of dimensions as: 

since = ~PAo is the state equation for a cosmological constant. 

Finally let us point out that the equations (1-10) may be expressed a, la Einstein by 
writing all extra terms due to the f{R) presence on the r.h.s. One can try to recover the 
standard form of the Einstein equations as follows 



C^. = R,.--g,.R^ -—^{T,, + T;lf) (1.18) 
where an effective energy-momentum tensor has been defined as 

(1.19) 



T^ff = - 



V,jR-\{f{R)-RfR)g^. 



This energy-momentum tensor does not necessarily obey the strong energy condition which 
holds in ordinary fiuids (dust matter, radiation, etc.) do. 
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1.4 Equivalence with Brans-Dicke theories 



Prom a classical field theory perspective, it is always possible to redefine the fields of a 

given theory in order to express the field equations in a more attractive way which would 
be easier either to handle or to solve. The price to pay is to introduce new auxiliary fields 
and even to perform either renormalizations or conformal transformations. 

It is widely assumed that two theories are dynamically equivalent if, under a suitable 
redefinition of either gravitational or matter fields, one can make the field equations to 
coincide. Nevertheless, some controversy has appeared in recent times especially when 
conformal transformations are used to redefine fields (see for instance [28] and [29] and 
references therein). 

As mentioned in Section 1.1 a possibility to construct alternative theories of gravity 
are the scalar-tensor theories which are based upon the introduction of an extra scalar 
field which modifies the gravitational sector. Those theories are still metric theories in the 
sense that the newly introduced fields do not couple to the fiuid contributions. 



The gravitational action for a general scalar-tensor theory in D dimensions is 



SsT 



d^a;v4 
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{d,<t>d^<t>)-U{<t>) 



2 2 

By choosing y(0) — (f)/K, uj{(j)) —uJo/{K,(j)) and U{(j)) — V{(t))/K, the action 
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d^xvT 
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ct>R-'^{d,ct>d^<j>)-V{cl>) 



:i.2o) 



(1.21) 



is obtained from (1.20). This is the action for the Brans-Dicke theories which is obviously 
a particular case of scalar-tensor theories. 

It can be shown that f{R) gravities within the metric formalism are nothing but a 
Brans-Dicke theory with Brans-Dicke parameter a;o = 0. This fact is easily proven as 
follows: a new field x is introduced and for the sake of simplicity let us define 



F{R) = R + fiR). 



[1.22) 



Thus the action (1.3) can be seen to be equivalent to the action 



^ 2k 



F{x) + 



dFjx) 



(R-x) 



[1.23) 



since if a variation of (1-23) with respect to x is performed, the equation which is found 
reads: 

dV(x), 



dx' 



-{R-X) = 



:i.24) 
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and thus x = -R provided d^/(x)/dx^ 7^ . Therefore the original action (1.3) is recovered. 
Defining now the scalar field as = dF(x)/dx and introducing a potential V{(f)) as 
follows 

y(0)=x(0)0-F(x(0)) (1.25) 
the action (1.23) takes the form 

which is exactly the same as (1.21) if o'o = is imposed. 

By including the corresponding fluid sector given by an energy-momentum tensor T^,^ , 
the field equations derived from (1.26) are 

= -| - ^a^Vi-P) + \'D^u<l> (1.27) 
R = ^ (1.28) 



where the trace of (1.27) 



{D - 1)00 + ^V{4>) + = ->^T (1.29) 

gives the dynamics of (j) in terms of the matter content. 

Let us finally note that if fuR = d^f{R)/dR'^ vanishes, the equivalence between the 
two theories cannot be guaranteed as can be seen from equation (1.24). On the other hand, 
the resulting Brans-Dickc equivalent theory makes clear that f{R) gravity theories have 
just one more extra degree of freedom than standard EH gravity. The apparent absence of 
kinetic term in the action (1.26) must not be thought of as the absence of dynamics in 
since this scalar is dynamically related to the matter fields, as can be seen from expression 
(1.29). Thus , or equivalently f{R) , is indeed a dynamical degree of freedom. 
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1.5 Geometrical results 

In this section we present different geometrical results obtained from the modified Einstein 
equations which were obtained in Section 1.3. Particular interest will be devoted in Sub- 
section 1.5.1 to vacuum solutions. Then, the possibility of mimicking the usual GR results 
using f{R) functions will be addressed in Subsection 1.5.2. These results were originally 
presented in [26]. 



1.5.1 Vacuum solutions 

Let us consider the EH action (1.2) in D dimensions with nonvanishing cosmological 
constant. In this case the equations (1-15) can be studied in vacuum, i.e. T^i, vanishes for 
all its components and therefore 

R^,,-^Rg^, + ^^-J^ADg^^. = (1.30) 

whose solutions satisfy 

R^u = ^Dg^u ; R = dAd (1.31) 

which motivated our choice for A^ normalization in Section 1.3. Equations (1.31) provide 
the conditions to be accomplished by a metric g^^, to allow vacuum solution in this case. 
If now one considers the f{R) general case provided by the equations (1.10), one may 
wonder about the condition for the existence of constant curvature solutions, i?o from 
now on, in a vacuum scenario. Thus, the equations (1.10) may be simplified to become 

R^. (1 + fR)-\ 9^.. {R + f{R)) = 0. (1.32) 

Note that the term involving V^^fR in (1.10) has disappeared since it vanishes when 
constant curvature is assumed. Taking the trace in the previous equation we get 

2R (1 + f^)-D{R + f{R)) = 0. (1.33) 

If Rq is a root of the previous equation, an effective cosmological constant may be defined 
as A^'^ = Rq/D . Provided the condition 1 + f'{Ro) 7^ is satisfied, Rq fulfills: 

Rq + f{Ro) 

Let us illustrate this procedure considering a simple model: 

f{R)^^ + go (1.35) 
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which has been widely studied in the hterature (see for instance [30] where D — A and 
go — 0). Then the constant curvature solutions - for an arbitrary number of dimensions 
D - are 

-Dgo ± ^D^{gl - Ag,) + 16^^ 

= 2{D^) ^ ^ ^ 

which reduce for D — A to the expression 

Ro^-go±^gi-Sg,. (1.37) 



For the EH case in D — A with cosmological constant A = A4 , i.e. gi — and 
go — — 2A , the constant curvature solutions are both Rq = 4 A and Rq = and for the 
vanishing cosmological constant case, i.e. go — , Ro = ±^/—3gl is obtained. 

As a different approach, one can consider equation (1.33) as a differential equation for 
the f{R) function so that the corresponding solution would admit any curvature R value. 
The solution of (1.33) is just: 

f{R) = aR'^/^ - R (1.38) 
where a is an arbitrary constant. Thus the gravitational action (1.3) becomes 

which has solutions of constant curvature for arbitrary it! . The reason is that this action 
is scale invariant since the ratio a/ k is a dimensionless constant. 



1.5.2 Some EH solutions reproduced by f{R) theories 

Now we shall address the issue of finding some general criteria to mimic, by using general 
f{R) gravities, some solutions of the EH action not necessarily of constant scalar curvature 
and either with or without a cosmological constant term. 

Let the metric tensor g^^ be a solution of EH gravity with cosmological constant, i.e. 
such that the equations (1.15) are fulfilled. Then the same metric tensor g^^, will be a 
solution for (1.10) provided the following compatibility equation 

fn R^, - i g^, [f{R) + {D - 2)Ad] + V,Jn = (1.40) 

is fulfilled. Note that the fluid content comprised in T^,^ has been considered to be strictly 
the same as in (1.15). This allowed us to cancel this term out in order to obtain the 
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compatibility equation (1.40). In Section 2.6, a slight deviation of fluid contents between 
EH and f{R) approaches will be permitted. 

Some particularly interesting cases in which to apply this approach are the following: 

1. The simplest case is obviously vacuum, i.e. T^^^ = 0, with vanishing cosmological 
constant A^, = . Then the equations (1-15) become: 

R^, = ^RQf.^. (1.41) 

which imply Rq = and R^^ = . Consequently g^i, is also a solution of any f{R) 
gravity provided the following condition 

/(O) = (1.42) 

is accomplished as seen from (1.40). This is for instance the case if f{R) is analytical 
around R = and it can be written as follows: 

oo 

f{R) = J]/„it:^ (1.43) 

n=l 

2. If the cosmological constant is different from zero {Ad ^ ), but still Tf^^ = , the 
constant curvature results given in (1-31) arc again obtained. Then the compatibility 
equation (1.40) reduces to (1.33) with Rq — DAjj. In other words, g^^ is also a 
solution of the f{R) case provided 

fiDAo) = AD{2-D + 2fR{DAn)). (1.44) 

Notice also that in this situation, i.e. nonvanishing and vacuum, according to 
the result in (1.39) there would also be a solution for any Rq in the particular case 
f{R) = aR'^/^-R. 



3. If the considered case is A/) = and conformal matter (T = T^^ — 0), then the 
equations (1-15) would imply 

i?o = ; R^, ^ -kT^, (1.45) 

which will have a metric tensor g^j, as solution. Therefore, provided 

/(O) = ; /^j(0) = 0, (1.46) 

the same g^i, is also a solution of any f{R) gravity. This result could have particu- 
lar interest in cosmological calculations for ultrarelativistic matter (i.e. conformal) 
dominated universes. 
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4. Again in the conformal matter case with nonvanishing A^i , constant curvature 

Ro = DAd (1.47) 

is a solution for (1.15) for a given metric Qfj^i, which is also a solution of f{R) pro- 
vided again that the condition (1.44) is satisfied. 

5. Finally for the general case with no assumption about Ad nor about T^,^ , the 
metric tensor g^^,^ will be a solution for any f{R) gravity but for a modified energy 
momentum tensor T^,^ given by: 

T^. = T^u-^l^fRR;.u-l[f{R) + {D-2)AD]g,^. + V^jRy 

(1.48) 



1.6 Constraints on f{R) theories to ensure viability 

f{R) gravity models turn out to be severely constrained in order to provide consistent 
theories of gravity. In this section we review both cosmological and strictly gravitational 
conditions presented in [31]. Some relevant bibliography will also be provided. 

The usual four conditions that are required for a viable f{R) theory are: 

1. Jrr > for high curvatures [32]. This is the requirement for a classically stable high- 
curvature regime and for the existence of a matter dominated phase in the cosmological 
evolution. In the opposite case, an instability, referred to in the literature as the 'Dolgov- 
Kawasaki' or 'Ricci scalar' or 'matter' instability, would appear. Indeed, if fuR is smaller 
than zero, then the extra degree of freedom of the theory would behave as a ghost. This 
stability condition may also be recovered in studies of cosmological perturbations [33] and it 
can be given a simple physical interpretation as in [34] where if an effective D dimensional 
gravitational constant is defined as Gg// = G'd/(1 -|- /r) then 

It is easy to notice from the previous equation that if fjm < 0, Ge// would increase as 
R grows since R itself generates larger and larger curvature via equation (1.12). Such a 
mechanism would act to destabilize the theory with no stable ground state since if a small 
curvature starts growing it will do so without limit and the system would run away. If on 
the contrary > , a negative feedback mechanism operates to compensate the growth 
of R and consequently the runaway behaviour will not appear ^. 



^Note that in this analysis l + fn has been supposed to be positive (i.e. Ge// > ) as will be required 
from the second condition below to ensure viability. 



18 



CHAPTER 1. INTRODUCTION TO MODIFIED GRAVITY THEORIES 



2- 1 + > for all Ricci scalar curvature values. This condition ensures the effective 
Newton's constant to be positive at all times as can be seen from equation (1.18) and 
the graviton energy to be positive. This condition will also be proven in Chapter 4 to be 
required to recover standard thermodynamics of Schwarzschild-anti-de Sitter BHs in f{R) 
theories. 

3- /i? < ensures ordinary GR behaviour is recovered at early times. Together with 
the condition Jrr > , it implies that /r should be negative and a monotonically growing 
function of R in the range — 1 < /r < . 

4. <^ 1 at recent epochs. This is imposed by local gravity tests [33], although it is 
still not clear what is the actual limit on this parameter and some controversy still remains 
about the required value [20, 35]. This condition also imphes that the cosmological 
evolution at late times resembles that of ACDM . In any case, this constraint is not 
required if one is only interested in building models for cosmic acceleration. 

Let us summarize this section by saying that viable f{R) models can be constructed 
to be compatible with local gravity tests and other cosmological constraints [36] . 

1.7 Brane- world theories 

As mentioned in the Motivation section, many of the SM extensions try to solve open issues 
in modern physics such as the hierarchy problem. Some approaches try to answer those 
questions by introducing extra spatial dimensions, where the number of dimensions of the 
total space (bulk space) D = 4: + S . Those attempts were first proposed independently by 
Kaluza [37] and Klein [38] and many proposals followed throughout the twentieth century 
[39, 40]. 

First proposals for large extra dimensions were provided in [40]: in this model, the SM 
matter is confined in a spatial 3-dimensional manifold and the brane itself is considered not 
to be a gravitational source. Hence the background metric is assumed to be Minkowskian. 
Gravitational fields are the only fields able to propagate through the whole bulk space. 
Therefore gravity also propagates in the extra d dimensions which are for simplicity often 
compactified in a toroidal shape whereby all extra dimensions acquire a radius Rb ■ 

One of the most important consequences of this hypothesis is the relation between the 
fundamental gravitational scale in D dimensions Md and the Planck scale Mp which is 
not a fundamental constant any more but the effective gravitational constant in the theory 
reduced to 4 dimensions. In fact one may write 

Mj> = VsMl+^ (1.50) 

where Vs is the compactified volume in 5 dimensions [40], for instance in the toroidal 
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case Vs — {^t^RbY ■ The expression (1.50) allows to reduce the fundamental gravitational 
scale to the electroweak scale, M/j ~ TeV, if extra dimensions are large enough. For 
instance, compactification scales of i?^^ ~ 10^^ eV to 10 MeV provide this effect for 
extra dimensions 5 ~ 2 and 7 respectively. By reducing the fundamental scale to Md , 
gravitational effects may be detectable in experiments involving energies of this order [40] 
as will be explained in the Section 1.9. 

1.8 Excitations in brane worlds: branons 

Since no relativistic object may be considered as rigid in relativistic theories, the 3-brane, 
when embedded in the total space-time, may present fluctuations. These fluctuations were 
originally studied in [41]. In the extra dimensions of the BW models, such fluctuations are 
usually referred to as branons. They give rise to new states whose low energy dynamics has 
been widely studied [42, 43, 44]. A vast bibliography can be found dealing with branons 
[45, 46], their predicted detection in future colliders experiments [47] and the explanation 
that they may provide for the origin of dark matter [25, 48]. 

Let us consider a D dimensional bulk space Md wherein the brane lies embedded 
and that for simplicity we shall assume to be factorized in the form M.^ = M.i x B where 
is a 4-dimensional space-time and B is a. 5 -dimensional compactified manifold. The 
brane is therefore assumed to he on the M.^^ space-time manifold. As already mentioned, 
the gravitational contribution of the brane itself will not be considered. 

Let us denote the coordinates over the manifold A^d as {x^, |/™} with /i, = 0, 1,2,3 
and m = 1, 2, 5 and the ansatz for the total space M.d bulk metric will be 



with signature — ; — , — ) . 

In the absence of the 3-brane, this metric possesses an isometry group that is assumed 
to be of the form G{M.d) — G{M.i) x G{B) . The presence of the brane spontaneously 
breaks the symmetry to some subgroup G(A^4) x H with H C G{B) some subgroup of 
G{B). Therefore the quotient space K = G{Md)/{G{M4) x H) = G{B)/H may be 



The position of the brane can be parameterized as Y^{x) = {x'^,Y"^{x)} where the 
first four coordinates of the total space have been chosen to be the space-time coordinates 
corresponding to the brane {x^} . Let us assume that the brane is located at a point on 
B , i.e., Yq = {Y'^{x)} corresponds to the fundamental state of the brane. In this case 
its induced metric in the ground state is just g^i, = Q/^,^ = G/^^, . However, when brane 




(1.51) 



defined. 
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Figure 1.1: Brane with trivial topology in Ais = M.2 x as originally presented in 
reference [49]. The fundamental brane state is plotted on the left whereas on the right side 
an excited state is presented. 

excitations (branons) are present, the induced metric becomes 

g^.. = d,Y'''d,Y''GMN = g,. - d.Y'^d.Y^g'^^. (1.52) 

This situation may be illustrated by the simple Figure 1.1 where a 1-brane (string) is 
represented within a total space with two spatial coordinates A^3 = A^2 x . 

Since the brane creation mechanism is in principle unknown, or at least out of the scope 
of the present section, let us assume that the brane dynamics is described by an effective 
action, so we are allowed to consider for this action the most general expression which 
is invariant under brane coordinates reparametrizations. Therefore, it is very common 
to perform an expansion in derivatives of the induced metric given by equation (1.52) to 
describe the brane dynamics. Then, the first order of this effective action would describe 
the brane dynamics at low energies and it is usually referred to as the Dirac-Nambu-Goto 
( NO ) action: 

Sng = - f j d^xVH (1-53) 

where a constant / with energy units appears, which may be identified with the brane 
tension r = and d^a;^/!^ is the brane volume element. As was mentioned above, the 
presence of the brane will break any existing isometry of B except those which leave the 
point Yq on B invariant. In other words, the group G{B) is spontaneously broken to 
H(Yq) denoting the Yq isotropy group. 

The brane excitations with respect to the broken Killing fields in B correspond to the 
zero modes and they are parameterized by the branon fields 7r"(x) , a = l,...,k where 
k = dim{G{B)) — dim(if ) . These fields 7r"(a;) may be interpreted like the corresponding 
coordinates in the quotient manifold K = G{B)/H . 
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In particular, for a fundamental state independent of the position Yq of the brane in 
the B space, the action of an element of G{B) over Yq will take Yq to the point on B 
with coordinates 

Y^{x) = y-(yo, TT^ix)) = V + C(^o)7r"(x) + 0{7r') (1.54) 

where the branons fields normalization is performed through k = 87r/Mp. At this stage 
it is important to stress that coordinates for the transformed point given by (1-54) only 
depend on 7r"(a;) , i.e., on the corresponding transformation parameters of the broken 
generators. 

If B is considered to be an homogeneous space, the isotropy group does not depend 
on the particular chosen point where the brane lies, i.e. -^^(^0) = H . In this case B is 
homeomorphic to the coset K — G{B)/H which is the space of the Goldstone bosons 
associated to the spontaneous isometry breaking - transverse translations - produced by 
the presence of the brane. Thus the transverse translations of the brane - branons - can 
be considered as Goldstone bosons on the coset K and the branon fields can be defined as 
coordinates tt" on K , which are chosen to be proportional to B coordinates, since the 
number of Goldstone bosons is equal to dim{B) , as: 

= ^Ci^'" (1-55) 
Kb 

where 

V = fRB (1.56) 

is the typical size of the coset K and Rb is the typical size, in length units, of the 
compactified space B . 

Therefore, according to the previous assumption (1.55), it is obvious that 



Qnu ^ Qtiv--^ Kp{'K)d^T:'^ d^T^^ (1.58) 



and the induced metric on the brane (1.52) is rewritten in terms of the branon fields tt as 

1 

where is the K metric which is easily obtained from the B metric 

BY"^ BY"- 

KM= f'~9'mniY{.))^^ (1.59) 

as explained in [45]. In more complicated cases in which translational isometrics in the 
bulk space are not only spontaneously but also explicitly broken, the metric g^^, could also 
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be a function of the extra dimension coordinates {y'^} . Then it is possible to show that 
branons may become massive. In fact in [25] these massive branons were shown to behave 
as WIMPs and thus form natural candidates for dark matter in this kind of scenario. 

Therefore for small brane excitations in a background metric , the effective action 
(1.53) can be expressed as a derivative expansion as follows: 

Seffin] = Si%[n] + 5g),[7r] + Si%[n] + ... (1.60) 
where the corresponding zeroth order is 

JM4 

Note that •S'^jjfTr] and S'gjj[7r] hold for contributions to the effective action containing 
two and four derivatives of the branon fields respectively. Let us finish this digression by 
remarking that the term S^^j:[7r] , with two field derivatives, is nothing but the non-linear 
sigma model action associated to the coset space K . 



1.9 Brane-skyrmions 

Apart from branons, the brane may support other states due to the nontrivial homotopies 
of the coset space K such as strings, monopoles or skyrmions. This fact appears due to 
the possibility of wrapping around the extra dimension space B giving rise to nontrivial 
topological configurations as was studied in detail in the reference [49] . 

In fact, texture-like configurations, called brane-skyrmions, arise when the third homo- 
topy group of K is nontrivial. In particular, for 

7T,{B) = n;{K) = Z, (1.62) 

the third homotopy group will be the minimal one supporting the existence of those non- 
trivial configurations^. Those brane-skyrmions can be nicely understood in geometrical 
terms as some kind of holes [50] in the brane which make it possible to pass through them 
along the B space. This is because in the core of the topological defect the symmetry 
is reestablished. In particular, in the case we are interested in, the broken symmetry is 
basically the translational symmetry along the extra- dimensions. 

In order to simplify the calculations, we shall consider an homogeneous compactified 
manifold B and the coset space K homeomorphic to SU (2) and equivalently to . 
Then 



Bc^K SU{2) ~ S^. 
^Note that 773 (B) = 773 (if) if B is an homogeneous space. 



(1.63) 
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Figure 1.2: Brane-skyrmion configuration with nw = 1 and nonvanishing size in A^3 = 
A^2 X as originally presented in [49]. 



This fact allows to hold a third homotopy group Z . 

Furthermore, let us introduce spherical coordinates on both spaces, A^4 and K as 
follows: in A^4 we denote the coordinates {t,r,9,ip} with G [0,27r), 9 G [0, vr] and 
r G [0, oo) . On the coset manifold K , the spherical coordinates are denoted {xk, ^K) <Pk} 
with (pK ^ [0, 27r) , 9^ G [0, vr] and xk ^ [0, vr] . Notice that such coordinates cover the 
whole spherical manifolds and relate to the physical branon fields (local normal geodesic 
coordinates on K ) by: 

TTi = V sinxK s^^^K cos (pK, 

712 = V sin Xk sin 9k sincpK, (1.64) 
TTa = f sin Xk cos 9^ ■ 

The coset metric in spherical coordinates is written as 

h^p= I v^sin\xK) I . (1.65) 

f2sin2(xx)sin2(^^) 



1.9.1 Static brane-skyrmions 



For static configurations, it can be proven that the mass for the brane-skyrmion may be 
obtained directly as: 

M[it] = - [ d'xCeff = f [ d'x,/\i\ (1.66) 

JMs JMz 
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where the effective Lagrangian comes from the expression (1.53). In general this expression 
is divergent due to the contribution of the zeroth order term when ^/\g\ is expanded in 
branon fields derivatives, refiecting the fact that the brane is an infinite object with finite 
tension. To prevent that, we substract this term in order to get a brane-skyrmion finite 
mass, i.e.: 



As was shown in the previous section, the vr" fields are mappings from the AI4 manifold to 
the coset manifold K . For static, i.e. time independent, field configurations, these could be 
understood as mappings from the corresponding spatial 3-dimensional hypersurface ( A^a ) 
to the coset space ( in the case we are studying). For finite energy configurations, fields 
should vanish at the spatial infinity and M.3 can be compactified to -S"^ . Therefore one 
may write tt" : 5^ — )> S*^ . Since the third homotopy group of is Z , the mappings 
can be classified by an integer number . Thus, branons can be identified with the 
topologically trivial configurations nw = , whereas those configurations with nw ^ 
will be denoted as brane-skyrmions. 

Consequently, for static skyrmions this mapping may be implemented in the following 
way: 



with the boundary conditions F{0) — F(oo) = n^yvr for a winding number nw 7^ . 

In this case, M^fvr] may be written as a F{r) functional and the correct mass for this 
kind of skyrmions is obtained by minimizing Ms[F] in the space function with adequate 
boundary conditions. 

From expression (1.67), it can be proven that the skyrmion is point-like, stable and its 
mass becomes: 




(1.67) 



(pK ^ (p ; Ok ^ 9 XK ^ F{r) 



(1.68) 



Ms = 27rV^i?| 



(1.69) 
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1.10 Gravitational signatures at the LHC 

As already mentioned, it can be seen from the equation (1.50), the fundamental gravita- 
tional scale could be as low as the electroweak scale if extra dimensions are large enough. 
Since the LHC is operating at a center of mass energy of ^/s = 14 TeV , if the fundamental 
scale of gravitation is Md ~ TeV , both the production and decay of Schwarzschild mini 
BHs at high ratio becomes possible [51]. 

These BHs once produced would decay into SM particles with a clean signature and 
a low background. Several features of such objects could then be extracted from expe- 
rimental data: for instance BH masses Mbh niay be determined very precisely due to 
the absence of missing energy and their temperature could be extracted from the energy 
spectrum of the products. Thus the correlation between these two quantities may provide 
relevant information able to determine the number of extra dimensions, and therefore the 
fundamental scale of gravity. On the other hand the Hawking evaporation law could be 
tested experimentally. 

The total cross section when two partons collide at the LHC with an impact parameter 
less than the Schwarzschild radius Rs is of order 

2 

D-3 

(1.70) 



and it does not contain small coupling constants. If Mn ~ TeV the cross section is of order 
TeV~ ~ 400 pb and therefore BHs will be produced copiously. The total production cross 
section ranges from 0.5 nb for Md = 2 TeV , D = 11 to 120 fb for Md = 6 TeV D = 7 . 
For Md ~ ITeV, the LHC - with a peak of luminosity of 30fb~^/year - will produce 
10^ BH/year. 

Experimental signatures rely on two qualitative properties: on the one hand, the 
absence of small couplings as seen from expression (1.70) and on the other hand, the flavor 
independence nature of BHs decays as will be explained in the following paragraph. Note 
that when Mbh approaches Md , some stringy corrections to the previous assumptions 
may arise but semiclassical arguments remain valid as long as Mbh ^ Md ■ 

Once the BHs have been produced they decay following a process governed by their 
Hawking temperature Th ~ 1/Rs with an associated wavelength A = 2tt/Th larger than 
the BH size and therefore BHs would emit, in a first approximation, as point radiators 
mostly in the s-waves. This indicates that BHs decay equally to particles on the brane and 
in the bulk since the decay is only sensitive to the radial coordinate. If the approxima- 
tion Mbh ^ Th is made, the average multiphcity of particles {N) produced in the BH 



(7(Mbh) ^''Rl=Ml 
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evaporation is given by: 



(N) 



D-3\Md 



2Vi /Mbh 




D-3 



(1.71) 



Since the decay is thermal, it does not discriminate between particle species (of the same 
mass and spin) and therefore BHs decay, roughly speaking, with the same probability to 
all SM particles. The signal of hard primary leptons and hard photons is quite clean with 
a negligible background since the production of SM leptons or photons occur at much 
smaller rate than BH production [51]. 

The way to determine M^u and Th deals with the study of decay products and the 
fits of the energy spectrum of those products to the Planck formula respectively. Once 
those two quantities are determined, they could provide some evidence of the Hawking 
radiation and of the fact that the observed events indeed come from BH evaporation and 
not from any other mechanism. 

The relation between those two quantities, Mbh and Th obtained independently, may 
shed light about the dimensionality of the space since it can be proved that 



where the constant does not depend on Mbh • Therefore the previous equation provides a 
direct method to determine the dimensionality D of the space as the slope of this relation. 

The experimental signatures outlined above allow us to state that if the fundamental 
scale of gravitation is of order TeV, as suggested in BW scenarios, some important physical 
consequences may appear. In fact, colliders study of BHs - eventually produced at a high 
rate in accelerators such as LHC - could help reveahng the main features of physics in the 
vicinity of the electroweak scale or even determining the total number of dimensions of the 
space-time. 



logiTn) 



1 



log(MBH) + constant 



(1.72) 
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Chapter 2 

Dark energy in f(R) theories 



2.1 Introduction 

As was commented in Chapter 1, when the modified Einstein equations were rewritten 
d la Einstein, the presence of a function f{R) in the gravitational sector modifying the 

usual EH Lagrangian may be understood as the introduction of an effective fluid which is 
not restricted to hold the usual energy conditions. Therefore f{R) functions may be used 
to explain the present cosmological acceleration. Historically, some f{R) models were 
proposed to modify GR at short scales, i.e., high energies trying to explain infiation, as for 
instance f{R) oc R^ , but no interest was paid in those models to provide a mechanism to 
cause late time acceleration. First attempts to induce cosmological acceleration considered 
f{R) (xl/R but those models turned out to be in conflict with solar system tests [52] and 
even to be unstable when matter is introduced [53]. 

Before studying these issues, let us mention that f{R) models, apart from satisfying 
those gravitational and cosmological conditions given in Section 1.6, should verify some 
extra conditions of cosmological viability. For instance, they have to include a background 
evolution with Big Bang nucleosynthesis (BBN) and both radiation and matter dominated 
cosmological eras. This fact will be explicitly studied in this chapter in Section 2.6. On the 
other hand, they must provide cosmological perturbations compatible with cosmological 
constraints from CMB and large scale structures (LSS). This fact will be studied thoroughly 
in Chapter 3. 

The present chapter is organized as follows: in Section 2.2 we shall revise the standard 
approach to describe the cosmological evolution in the ACDM model in a homogeneous, 
isotropic and spatially flat metric. In the following Section 2.3 we shall generalize the 
usual Einstein equations when f{R) gravity theories are present. Then we shall study in 
Section 2.4 the cosmological viability conditions for f{R) theories to hold a dust matter 
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dominated era followed by a late time acceleration and some interesting models which have 
been proposed to be viable will be provided in Section 2.5. Then, Section 2.6 will be the 
core of the chapter and it will be devoted to study if /(-R) theories are able to mimic 
standard ACDM evolution. These f{R) models will possess vacuum solutions with null 
scalar curvature what allows to recover some GR solutions usually considered. 

To finish this chapter, we shall study in Section 2.7 how the modification of the gravi- 
tational sector by f{R) models may mimic the infiuence of perfect fiuids (parameterized 
by a constant equation of state) in the cosmological evolution without any presence of 
such a fluid in the fiuid content. The chapter will finish with Section 2.8 by drawing some 
attention over the main obtained conclusions. 

The results presented in this chapter were originally published in [54]. 



2.2 Standard Einstein equations in a FLRW universe 

Since the leitmotiv of this chapter is to study cosmological solutions, our universe, which is 
assumed to be isotropic and homogeneous at large enough scales for fundamental observers, 
may be represented with a D = 4 dimensional Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metric 

/ dr^ 



ds' = dt' - a\t) + r'dn'^j (2.1) 

expressed in cosmic time t and where a{t) is usually referred to as the scale factor. 
Alternatively, this metric may be expressed in conformal time r , defined by the relation 
dt = a(r)dr and thus this metric becomes 

ds^ = a\T) (^dr^ - - r^dnl^ (2.2) 

In this metric, the Hubble parameter may be defined in either cosmic or conformal time 

as 

^ da{t)/dt ^ a Mrydr^a^ 
^ ' a{t) a a{T) a{T) ^ ' 

respectively and the identity aH = is straightforwardly inferred. 

For the values of the parameter k smaller, equal or bigger than zero, the universe is 
spatially hyperbolic, fiat or spherical respectively. In the following calculations we shall be 
considering /c = . This choice is justified according to WMAP data [11] where the results 
obtained for ACDM model are: = -k/Hl, with Hq = H{ttoday) = 100/ikms~^Mpc~^ 
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and h = 0.699 ± 0.018 and -0.0133 < < 0.0084 (95% CL). Therefore terms related 
with k will be subdominant in either Priedmann's or generalized Priedmann equations to 
be presented in following the section. 

Considering the previously introduced metric (2.1) and perfect fluids given by (1.7) 
for the present fluids, the only two independent Einstein equations for D = 4 are the 
Friedmann and the acceleration equations respectively, which may be written in cosmic 
time as: 



- = — J]p« + 3P« 2.5 

a 

where G = G4 is the gravitational constant in four dimensions and the summation over 

subindcx a: holds for the present fluids contributions (baryons, radiation, dark matter, 
DE, etc.). If positive cosmological acceleration is required, i.e. a > , the condition to be 
accomplished from expression (2.5) would be X]a(Pa + 3Pa) < 0. This condition would 
require that if only a perfect fluid is present, its state equation would satisfy the condition 
uJa < —1/3 which is not the case for standard fluids, such as for instance dust matter and 
radiation. On the contrary a cosmological constant does provide positive acceleration in 
equation (2.5) since its state equation satisfies uj\ — —1 . 

For this metric, the energy-momentum conservation equations lead, in cosmic and 
conformal time respectively, to the following equations: 

Pq, + 3(1 +UJa)H Pa = 

p'a + 3{l + iUa)npa ^ (2.6) 

which hold separately for each fluid whose state equations are assumed to be Pa — oUaPa ■ 
Previous equation is integrated to give: 

Pa{t)=Pa{to)(^X^'^^"^ (2.7) 



a{t)J 

where to is an arbitrary time and the corresponding scale factor is a(to) • 

By using the definition given in (1.5) and (1.6), the Ricci scalar curvature for FLRW 
spatially flat metric is written in terms of the scale factor and its derivatives as follows 



R = 6 



2 

a \ a 



a 



^(n' + n^). (2.8) 
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Let us finish this section by rewriting the previous Priedmann equation (2.4). To do 
so, let us divide that equation by Hq = H'^ito) and consider as present fluids dust matter 
{um = 0), radiation (wRad = 1/3) and cosmological constant A (a;A = — 1). Thus, 
taking into account (2.7) for each present fluid, we get: 



QMa(t)-' + QRada(^)~^ + ^^A (2.9) 



where we have used the notation: 

_ ST^Gpuito) _ 87rGpRad(to) A 

and the normalization of the scale factor a(to) = 1 has been used. Note that if expression 
(2.9) is evaluated at t — to, then Qm + ^Rad + = 1 ■ 



2.3 Modified Einstein equations in a FLRW universe 

Inserting the metric (2.1) for D = 4 in the equations (1.10) and assuming also energy- 
momentum tensor as given in (1.7) for a fluid with energy density po and pressure Pq , 
the only independent modifled Einstein equations are 

3{l + fR)--hR + f{R))-3-RfRR = -SnGpo (2.11) 
a 2 a 

(l + fj,)(H + 3H')-^{R + f{R))-^^^(^a'RfnR) = SnG P^ (2.12) 
and in conformal time r , these equations are given by 

-^{^ + fR)-l,{R + fm--YfR = -S^^Po (2.13) 

Qj ^ Qj 



\{n' + 2U^){l + fn)-UR + fm-\WR + fR)-^^GP,. (2.14) 

Remind that dot denotes here derivative with respect to time t whereas r derivative was 
denoted with prime. A very useful equation to use in the following calculations is the 
combination (2.14) minus (2.13) which becomes 

2(1 + fR){-%' + + 2%fj, - = 87rG'(po + P^)a\ (2.15) 

At this stage we should note that, for instance, according to equations (2.11), (2.12) 
together with (2.8), it is clear that modified Einstein equations are not second order in 
derivatives any more, but at least third order in the scale factor derivatives provided 
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2.4 Cosmo logical viability for f{R) dark energy mo- 
dels 

In this section we shall revise the conditions that a model for DE given by a f{R) theory 
must fulfill in order to be cosmologically viable: i.e., any viable f{R) model should have 
a matter dominated phase long enough to provide the adequate cosmological evolution 
prior to a late time acceleration phase. As a matter of fact, equations (2.11) and (2.12) 
with dust matter as the unique present fluid, can be rewritten in the form of a system of 
autonomous equations [55]. In that reference two variables, m and r are introduced as 
follows: 

^ RfRR . ^ R{1 + fn) ^2 16) 

i + ' R + fiR)- ^ ^ ^ 

Both the dynamics and stability of that autonomous system are determined by six critical 
points Pi,. ..,6 ~ according to the notation in [55] - that appear in the system resolution. 



2.4.1 Critical points and stability 

According to the results presented in [55] , both points P5 and Pq satisfy 

m(r) = -r - 1. (2.17) 

If in the previous equation m is assumed to be constant, the condition (2.17) holds straight- 
forwardly from two other equations in the autonomous system. In this case the points P2,...6 
always exist while Pi and P4 are present for values m — 1 and m — —1 respectively. 
The critical points P5 and Pq which give the exact matter era evolution, i.e., a{t) oc t^^^ , 
exist only for m = ( P5 ) or for m = — (5 =b v^73)/12 ( Pq ) but the latter corresponds to 
a vanishing matter density and obviously it does not give a standard matter era. 

If on the contrary m is not assumed to be constant, the number of solutions depends 
on the particular f{R) choice, but only Pi, 5,6 can be accelerated and only P5 might give 
rise to matter era. This last situation would require m ~ to resemble the standard 
matter era evolution. Summarizing the result in this case, only trajectories passing near 
P5 with 171)^0 at r^— 1 and landing on an accelerated attractor would give a viable 
cosmological evolution. 



2.4.2 Classification of f{R) models 

By studying all possible trajectories in the m and r variables of the already mentioned 
autonomous system, it can be shown [55] that a classification of f{R) models can be 
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based entirely upon geometrical properties of the curve m(r) . These two variables allow 
to classify the f{R) models in four different classes: I, II, III and IV depending on the 
existence of a standard matter epoch and a final accelerated expansion as follows: on the 
one hand, the viable matter dominated epoch requires 

din 

] — > -1 (2.18) 

dr 

at r ^ —1 . On the other hand, the late time acceleration epoch requires to fulfill one of 
the following two conditions: a de Sitter acceleration follows the matter epoch if and only 
if 

1. < m(r) < 1 at r = -2 (2.19) 

whereas a non-phantom accelerated attractor follows the matter dominated epoch if and 
only if 

2. m = -r - 1 ; ^^--^ < m < 1 ; < -1. (2.20) 

2 - ' dr ^ ^ 

For instance, according to the previous requirements over m and r variables, models 
of the type f{R) = aR^^^ — R and f{R) = aR^"' do not satisfy these conditions for any 
n > and n < —1 and are consequently cosmologically nonviable. 

The main features of each class of models are: 



Class I: this class covers all the cases for which the curve m(r) does not connect 
the accelerated attractor with the standard matter point (r, m) = (—1,0) either 
because m(r) does not pass near that matter point, i.e., m{r —1) 7^ , or be- 
cause the branch of m(r) that accelerates is not connected with the standard matter 
point. Moreover, instead of having a standard matter phase given by a scale factor 
a{t) oc t^/^ , these f{R) models possess a peculiar scale factor behaviour a{t) oc t^^^ 
before accelerating epoch and are therefore unsuitable models. 

Class II: for these f{R) models the m(r) curve does connect the upper vicinity 
(m > 0) of (r, rri) = (— 1, 0) with a critical point able to provide acceleration. There- 
fore models here have a matter epoch and are asymptotically equivalent (hardly 
distinguishable) to ACDM model {oOeff = ~1 — 2H /3H^ = —1), i.e., they are 
asymptotically de Sitter and observationally acceptable. These models satisfy both 
equations (2.18) and (2.19). 

Class III: these f{R) models may possess an approximated matter era but as a 
transient state followed by a final and strongly phantom attractor at late-time. This 
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f{R) models 


m{r) 


Class I 


Class II 


Class III 


-R + aR-"^ 


— 1 — n 


n > -0.713 




-1 < n < -0.713 


aR-"" 


n(l+r) 


n > 


n G (-1,0), a < 




-R + RP [log {aR)f 


qr 




p= 1, g > 




-R + RPexp qR 


-r + 2 








-R + RPex.-p{q/R) 


p+r{2+r) 
r 




p = l 





Table 2.1: Classification of some f{R) DE models presented in [55]. None of these models belongs 
to Class IV. Models that belong to Class II for the provided parameter intervals, at least satisfy the 
conditions to have a matter era followed by a de Sitter attractor. 



is due to the fact that the m{r) curve intersects the critical line m{r) — —r — 1 at 
—1/2 < m < 0. The approximated matter era is a very fast transitient phase and 
only a narrow range of initial conditions may allow it. Since matter era is practically 
unstable, these models are generally ruled out by the observations. 

Class IV: for models of this class the connection between the upper vicinity of the 
point (r, m) = (—1,0) to the region located on the critical line m(r) = — r — 1 
is possible. Therefore these models are observationally acceptable: they possess an 
approximate standard matter epoch followed by a non-phantom acceleration with an 
effective equation of state cjg// = — 1 — 2H/3H^ > — 1 , thus these models posses a 
standard DE behaviour. These models satisfy both equations (2.18) and (2.20). 

Classes II and IV have therefore some chance to be cosmologically viable but the 
basin of the attractor has to be determined to provide acceptable trajectories according to 
the already mentioned analysis fully performed in [55]. In Table 2.1 the previous analysis 
have been applied to some f{R) models usually considered in the literature. 



2.5 Some cosmologically viable f{R) models 

In this section we provide three f{R) models already presented in the literature which 
claim to be cosmologically viable. 

a)f{R)^XRo 

This model was originally considered in reference [56] with n, A > and Rq of the 
order of the presently observed effective cosmological constant. Then /(O) = and the 
cosmological constant is claimed to disappear in fiat space-time but /flij(O) is negative 
and therefore, according to condition 1 in Section 1.6, flat space-time would be unstable. 
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For scalar curvatures R ^ Rq , f{R) tends to XRq and the model would behave as 
the EH case with an effective cosmological constant. On the other hand, dc Sitter space- 
time with curvature i?i > is also a vacuum solution provided -Ri(/r(-Ri) — 1) = 2/(i?i) 
according to equation (1.33) for D = 4 and R = Ri . Thus, this case would present an 
effective cosmological constant A(i?i) = i?i/4. 

For this model it can be proved that conditions 1 and 2 of Section 1.6 are satisfied 
in the curvature interval oo) if they are accomplished at R = Ri . Then these two 
conditions hold if and only if 



1 + 



Ri 

Rn 



n+2 2 4 

>l + (n + 2)(^] +{n + l){2n + l)(^] . (2.21) 



Note that condition 3 in that section is straightforwardly satisfied if as considered, pa- 
rameters A and n are positive. On the other hand, this model also satisfies the required 
conditions to provide a matter dominated era at R^ Rq and does not possess the already 
mentioned Dolgov-Kawasaki instability. The remaining condition 4 is easily accomplished 
if Rq is considered much smaller than R at recent epochs. 

A simple choice of parameters A , n and Rq shows that this model obeys the conditions 
(2.18) and (2.19) and therefore according to the analysis in the previous section it possess 
a matter dominated epoch and a de Sitter late time acceleration. 

b) fiR)^-ami{^r[l + P{^)r' 

This model was first proposed in reference [33] in order to mimic ACDM evolution in 
the high-redshift regime and to accelerate at low rcdshift with an expansion history close 
to ACDM model. In this model, the parameter n was considered to be positive and for 
convenience the mass scale a was given by 

„^5!L|A. = (8315Mpc)-(^) (2.22) 

where Pq is the average density today. This model does not have a bare cosmological 
constant since /(O) = and its parameters mi , (5 and n may be rewritten as [57] 

777,1 1 — 0,M 

-i ^ 6—-^ (2.23) 



mi fniRo) ( 12 



?2 



n \Vt 



I'M 



n+l 



9 (2.24) 



where VIm is the effective matter energy density at the present time. Finally the constraint 
|/i?(-Ro)| < 0.1 was imposed and Rq is the scalar curvature today as would be obtained 
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from ACDM model, i.e. 




(2.25) 



Setting the values n = 1 and Qm = 0.3 , it was proven in reference [57] that this model 
belongs to the Class II presented in the previous section. 



This two-parameter f{R) model presented in [58], where parameters a and are 
positive, has claimed to be cosmologically viable and different from ACDM . In fact, it 
does satisfy both cosmological conditions (2.18) and (2.19) presented in previous Section 
2.4, provided a > 1 and regardless of the value of i?* . 

Conditions 1, 2 and 3 given in Section 1.6 are also satisfied but concerning the condition 

4 also given in that section, it will be shown in Section 3.4 that this f{R) theory docs not 
satisfy this condition, showing that this model is indeed distinct from ACDM . This fact 
and its consequences will be studied in Section 3.4. 



Now that the modified Einstein equations have been presented for FLRW metric, some 
interesting results at cosmological scales will be obtained in this section. The presented 
approach tries to mimic cosmological well-known OR results in different cosmological eras 
employing adequate f{R) functions. Reconstruction procedures of this kind have been 
widely studied in the literature [59, 60, 61] where by rewriting the involved equations in 
new variables and assuming a given cosmological solution, mainly in vacuum, the required 
f{R) gravity is obtained. 

In this section the addressed issue will be to find a f{R) gravity able to reproduce 
the current cosmic speed-up appearing in standard ACDM cosmology. This function is 
required to be analytical at R — and to have R — as a vacuum solution, therefore 
it will not contain any cosmological constant contribution. From a more formal point 
of view we are seeking for some f{R) gravity model having the same FLRW solution 
as the standard EH action with cosmological constant for nonrelativistic matter (dust 
matter, i.e. Pm = ). For the searched cosmological constant absence it is clear that the 
f{R) expansion at i? = must start at the R^ term to avoid, on the one hand, having 
cosmological constant and, on the other hand, redefining the gravitational constant. 



c) /(i?) = -ai?.log(l + f) 



2.6 f{R) with no cosmological constant 
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2.6.1 Cosmological evolution in ACDM model 

Let us solve Einstein's equations (1.15) for the standard EH action plus a cosmological 
constant A with dust matter in the energy-momentum tensor side. 

The most recent cosmological data quoted in reference [11] are compatible at late times 
with a cosmological model based on a spatially fiat FLRW metric like (2.1) together with 
Einstein's equations with a cosmological constant A 7^ and dust matter (including dark 
matter). In this case, the equations (1.15) will be vahd and matter content will be written 
in terms of a pressureless perfect fluid 

Tt= diag(pMo(t), 0, 0, 0). (2.26) 



Equation /i — u — (time-time component) in (1.15) becomes 



a„(t) J 3 3 

where pMo{t) in previous equation is given by expression (2.7) if uJa=M is fixed to 0. Thus 



PMo{t) = PMoito) (^1^) (2-28) 

where the subindex means that the standard EH equations (1-15) with a cosmological 
constant are being considered. This notation will be relevant later on when standard 
ACDM cosmology will be compared with the results coming from the action (1.3) for the 
f{R) function that we shall find in this section. 

Substituting the expression (2.28) in the equation (2.27), the solution ao{t) , using the 
notation introduced in Section 2.2 in this chapter, is found to be: 

ao{t) = (^^y\inh'/'(^-^Hot]. (2.29) 

On the other hand, by taking the trace of (1-15) in this case, i.e. cosmological constant 
and dust matter, it is found that 

Ro{t)-4A ^ 87rGpMo{t). (2.30) 



2.6.2 /{R) case with no cosmological constant 

Now let us consider the equations (1-10) but in the case where no cosmological constant 
is considered and the energy-momentum tensor for dust matter will be 

= diag(pM(i), 0, 0, 0). (2.31) 
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Then the equation (2.11) becomes 

3(1 + fR)--l:{R + fiR)) - ^-RfRR = -SnGpM (2.32) 

Qj ^ a 

where we have ehminated the subindex in the different quantities to avoid any confusion 
with the previous case presented in Subsection 2.6.1. As was already mentioned, it is clear 
that the solutions for equation (2.32) will strongly depend on the function f{R) : different 
choices for this function will lead to different evolutions of the universe for the same initial 
conditions. However, our approach to the problem will be to find a function f{R) so that 
the solution a{t) of the equation (2.32) will be exactly the same as the solution provided 
by the expression (2.29) that we obtained by using GR with nonvanishing cosmological 
constant and which seems to fit the present cosmological data. In other words, we want 
to find the f{R) model such that the solution for the equation (2.32) is exactly the scale 
factor (2.29), i.e.: 

a{t) = ao{t) (2.33) 

for the same initial (or present, i.e. t = to) conditions. If it were possible to find such 
a function f{R) then, it would be possible to avoid the necessity for introducing any 
cosmological constant just by considering a gravitational action such that given in the 
expression (1.3). In the following it will be shown that such a function happens to exist and 
its precise form will be provided. In order to do that one first notices that accomplishing 
the condition (2.33) after radiation-matter equality clearly implies 

R = R{t) = Ro{t) (2.34) 

and then R{t) and Ro{t) may be used indistinctly. On the other hand we shall write 
the matter density as the former matter density provided by expression (2.7) plus a new 
contribution, i.e. 

PM{t) = pMo{t) + Ap{t), (2.35) 

accounting for a slight variation with respect to the density provided in the Subsection 
2.6.1. Assuming that matter for arbitrary f{R) is still nonrelativistic in this cosmological 
era we have 

Ap{t) = Apito) l^^y (2.36) 

where according to the expression (2.7) particularized for the constraint (2.33) and consi- 
dering the relation given in (2.30) we can write 



/ao(to)y^ R{t) -4A 
\ao{t)j SttG pMoito) 



(2.37) 
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and then (2.36) becomes 

, , , R(t) - 4A , , 

Ap(t) = -rj^ (2.38) 

K 

where we have introduced the parameter 

PMo[to) 

SO that matter density (2.35) is rewritten as 

PMit]v) = {'^-v)pMo{t) (2.40) 

Finally the last term on the l.h.s. of equation (2.32) can be written in terms of the scalar 
curvature by differentiating expression (2.30) with respect to cosmic time t and using the 
conservation equation (2.6). Hence, we get 

3{R- 3A) {R - 4A) f^R + f- Ji? + 3A^ Jr - lf{R) - A 



-77(i?-4A) = (2.41) 

where the time dependence of R is implicit. This last equation can be considered as a 
second order linear differential equation for the function f{R) , so two initial conditions are 
needed to solve it: the natural choice that has been judged more convenient and physically 
meaningful is the following: 



1. Firstly, the absence of any cosmological constant in the gravitational action is re- 
quired, so that /(O) = . 

2. Secondly, the standard EH action behaviour should be recovered for low scalar cur- 
vatures without redefining the Newton constant, i.e. /r(0) = . 



Moreover, f{R) function is wanted to be analytical at the origin so that R — should 
be a solution for the field equations in vacuum. This is an extremely important requirement 
since it allows both Minkowski and Schwarzschild solutions to be vacuum solutions. 

With these initial conditions, the equation (2.41) can be solved by using standard 
methods. One particular solution is: 

UiR) = -r)R + 2ir)-l)A. (2.42) 

The homogeneous equation associated with (2.41) is a Gauss-type equation solved in terms 
of hypergeometric functions. The general solution of the homogeneous equation can be 
written as: 

MR) = A (KMR) + K-f-{R)) (2.43) 
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where 



-a± 



R 
A 



) 



-1 



2F1 a±, 1 + a± — c, 1 + a± — a=p; 



3 



(2.44) 



and the symbol 2F1 holds for hypergeometric functions and the constants 



a± = 



12 



{7±V73} ; c 



1/2 



(2.45) 



have been introduced. The rj -dependent constants K^^ and must be determined from 
the initial conditions given above. Numerically it is found that: 

= 0.6436 (-0.9058 77 + 0.0596) ; X_ = 0.6436 (-0.2423 77 + 3.4465). (2.46) 

The hypergeometric functions given in (2.44) are generally defined in the whole complex 
plane. However, a real gravitational action is wanted. In principle this requirement is very 
easy to achieve since the coefficients in the equation (2.41) and the constants K± are all 
real. Then it is obvious that the real part of the functions appearing in the expression 
(2.44) is a proper solution of the homogeneous equation associated with (2.41). Thus the 
function we are seeking can be written as: 



Nevertheless, the situation is more comphcated. The homogeneous equation has three 
regular singular points at Ri = 3A , R2 = 4A and R^ — oc . This results in the 
solution fh{R) having two branch points -Ri and R2 ■ More concretely there are two cuts 
along the real axis: one from minus infinity to Ri and another from R2 to infinity. Thus 
one must be quite careful when interpreting (2.47). Prom minus infinity to Ri there is 
only one Riemann sheet of fh{R) where /(O) and /ii(0) vanish and therefore this is the 
one that we have to use to define f{R) . From Ri to R2 the real part of fhiR) is well 
defined. Finally from R2 to infinity there is only one Riemann sheet producing a smooth 
behaviour of f{R) ■ To reach this sheet one must understand R in the above equation 
as R + ie . At the present moment we do not know if this analytical structure has any 
fundamental meaning or it is just an artefact of our construction. Much more important 
is the fact that the function R + fp{R) + fh{R) , which is the analytical extension of our 
Lagrangian, is analytical at R = , having at this point the local behaviour R + 0{R^) . 
Therefore our generalized gravitational Lagrangian R+f{R) does guarantee that i? = is 
a vacuum solution as can be seen from expression (1.12). At the same time, this Lagrangian 
reproduces the current evolution of the universe without any cosmological constant. 

Now that the f{R) function in (2.47) has been obtained, it is possible to check our 
result out by solving (2.11) in terms of a{t) for the f{R) given in (2.47). This is done 
by rewriting the equation (2.32) in terms of a{t) by using (2.8) and (2.7) together with 
the dust matter density in terms of 77 and a{t) as provided by expression (2.35). At this 



f{R) = U{R) + Re[hm- 



(2.47) 
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stage, the consistency of our results has been checked out by introducing the scale factor 
(2.29) and f{R) given in (2.47) in (2.11) and both sides of the equation turned to be equal 
for all 7] values. Thus it has been guaranteed that our gravitational action proportional 
to R + f{R) provides the same cosmic evolution -in the required cosmological eras - 
as the EH action with cosmological constant A in a dust matter universe. Therefore, 
our model does verify, in the same range of precision, all the experimental tests that the 
standard cosmological model does in the present era. Notice also that in principle this can 
be achieved for any value of r] , i.e. for any desired amount of matter. Nevertheless, some 
restrictions should be imposed on the parameter rj . For instance it is obvious that in a 
dust matter dominated universe Pm(^; ^) > implies t] < 1 . 

Much more stringent bounds can be set on parameter rj by demanding that this model 
works properly back in time up to Big Bang Nucleosynthesis (BBN) era. Observations in- 
dicate that the cosmological standard model fits correctly primordial light elements abun- 
dances during BBN , that means that the expansion rate H{t) cannot deviate from that 
of standard cosmology i?o(^) in more than 10% for the background evolution (see for 
instance [62] for further details). Therefore by the time of BBN, departure of our model 
from the standard cosmology must not be too large and the equation (2.11) should provide 
a similar behaviour to the one given by the standard Priedmann equation (2.27) where 
now the density will include both dust and radiation contributions. 

At BBN the DE contribution is neghgible compared with dust and radiation densities. 
The scalar curvature is of order 10~^^eV^ (with h — c — 1 ior these calculations) 
and by that time dust and radiation densities are of the order of lO^^eV^ and eV"^ 
respectively. Since ~ i?o we can rewrite the equation (2.11) as a modified Friedmann 
equation as follows 



As was commented above, to reproduce light elements abundances it is required that 
if^(t) = i/o(t)(l ±0.2) for curvatures of order i?BBN • This implies that the second factor 
on the r.h.s. of the previous expression (2.48) should be between 0.8 and 1.2 by that 
period. Thus in order to match our f{R) gravity model with the standard cosmology at 
the BBN times we need to tune r] to a value about 0.065 with a stringent fine tuning. 
Therefore the matter content of our model is not too different from the one in the standard 
cosmology and the difference is in fact smaller than experimental precision in [11]. 

Concerning the problem of viability for this particular f{R) model, if conditions 2, 
3 and 4 given in Section 1.6 are required to hold, the parameter rj has to be fixed to 
a fine tuned value rj ~ —1.4311 but for this value fun reverse its sign at high enough 
curvatures and therefore the condition 1 in that section is not accomplished. Therefore it 
may be stated that the f{R) model given by expression (2.47) should be considered as an 




(2.48) 
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effective model able to reproduce ACDM model cosmological expansion after radiation- 
matter equality but not as a consistent gravitational theory valid for all scales. 



2.7 Effective fluid description of f{R) gravities 

In this section we shall be interested in finding those f{R) functions such that the corres- 
ponding modified Einstein equations in vacuum exactly reproduce the cosmological evolu- 
tion of EH gravity with a given perfect fiuid, i.e., the introduced modifications of the EH 
action through a function f{R) will play the role of the fiuid source. 

Let us thus consider such a perfect fiuid 'f obeying the following barotropic state 
equation 

Pf = Ufpf (2.49) 

with constant Uf and whose density pf scales according to the conservation equation (2.7) 
with the scale factor a = a{t) as 

p,{a) = p,{a{to)) (^^j = pMto))x'+^^ (2.50) 

where pf{a{to)) is the value of the fiuid density for a given value of the scale factor a{to) = 
1 . For the sake of simplicity a new variable x has been introduced in the previous 
expression (2.50) defined as follows 

X ^ 1. (2.51) 

With this new variable x , if the only present (or at least the dominant one) fluid is the 
deflned above, the standard EH Friedmann and acceleration equations in cosmic time t 
are respectively 

o 

- = -^p,{l + 3uj,) ^ -l-H^{l + 3u,)x'+^^ (2.52) 

where as usual 87rG'pf(a(to)) = Si^g . 

Analogously to the procedure in the previous section of this chapter, we shall consider 
the modified Einstein equations where there will not be any fiuid contribution. The solution 
for these equations is wanted to be the same scale factor as the one which is the solution 
of equations (2.52). In other words, the presence of the fiuid in the EH equations wants to 
be replaced by the contribution of some f{R) function in the modified Einstein equations. 
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To do so, note that the scalar curvature according its expression (2.8) for spatially flat 
FLRW metric may be rewritten as 



and 



R = 6 



dR 



a\ a 
- +- 
a J a 



871^(1 - 3a;f)pf 



R= — = -3H87rG{l + uJi){l-3uj{)p{{a{to)). 



(2.53) 



(2.54) 



At this stage, let us introduce a dimensionless variable R = R/Hq. Hence with this 
notation in variables R and x we get from expression (2.53) that 



/ 1 \ 3(l+a)f) 

^ = 3(1 - 3a;f) f - j = 3(1 - 3uji)x^+^'. 



(2.55) 



Note that for the specific choices Uf = —1, 1/3 no relation between variables R and 
X may be straightforwardly established through (2.55). For those cases x has to be 
determined by solving the equations (2.52). 



Since x — —3 Hx the following equalities can be written down 



= X 



l+£Jf 



(2.56) 



R = -9H (1 - 3a;f)(l + ujf)x^+'^' 



(2.57) 



d R 

— = R^ ^ 3(1- 3uji)(l + uj{)x'^' 
dx 



(2.58) 



d / 1 
\R, 



-i-3(l - 3ujf){l + uj{)ujfx'^'-^ 



(2.59) 



d/ ^ d^J_ . ^2 
d^ dx r' ° 



di?2 



d'fix) 1 d/(a;) 3 



da;2 dx Rl 



^{l - 3ui){l + Ui)uix'^'-^ (2.60) 



—R — — 3 H^xRx Hn 



(2.61) 
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where f{R) = f{R)/H^. 

Therefore, by considering a fluid density given by expression (2.50), the modified Fried- 
mann equation given by (2.11) may be rewritten as a second order differential equation 
for f{R) . Analogously, this equation may be expressed as a differential equation in x 
variable for f{R{x)) = f{x) and it becomes 



dfjx) 1 
dx R^ 
= 3H^x^+''' 



.a ~ 



3-R^ - 9H^ [3(1 - 3a;f)(l + uji)ujfx'^'] } - -H^ f{x) + 9H 



x d'^f{x) 
Rx dx'^ 

(2.62) 



where expressions (2.52) and (2.58) have to be substituted in the previous equation. 
The general solution of (2.62) is then: 

f{x) = -3(1 - 2,uji)x^^'^' + c+x'^t + c.x'^t" (2.63) 

where 



± 



1 

12 



9 a;f + 7 ± (9^;^ + TSc^f + 73) 



1/2 



(2.64) 



Requiring 7^ —1, 1/3 to avoid possible indeterminacies, it is possible to rewrite (2.63) 
in terms of R as 



f{R) ^ -R + C+ 



R 

3(1 - 3cjf) 



+ c_ 



R 

3(1 - 3cc;f) 



l+Wf 



(2.65) 



2.7.1 Some examples 

Some interesting cases for the fluid content are the dust matter, radiation and cosmological 
constant fluids, , i.e., oom, Ra.d,A = 0,1/3,-1 and rji = ?7M,Rad,A respectively. For these 
three cases, the corresponding functions become 

fuix) = c+xU'+^^) +c_xU'-^) -3x 
fj,{x) = -12 + c_x-^/2 + c+ 
fKM{x) = c+x'/' + c^ (2.66) 



where constants c± are arbitrary integration constants that can be fixed if either boundary 
or initial conditions are imposed. 
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2.8 Conclusions 

In this chapter we have found the f{R) gravity which exactly reproduces the same evolu- 
tion of the universe, from BNN up to the present time, as standard ACDM model does, 
but without the introduction of any form of DE or cosmological constant. The gravitatio- 
nal Lagrangian R + f{R) is analytical at the origin and consequently R = is a vacuum 
solution for the field equations. Therefore Minkowski, Schwarzschild and other important 
= GR solutions, with A = , are also solutions for this f{R) gravity. This result 
was originally presented in [54]. 

The price that we have to pay for all those good properties is that our Lagrangian, 
considered as a function of R , has a very complicated analytical structure with cuts 
along the real axis from infinity to R = 3A and from R = 4A to infinity. Obviously 
the only reasonable interpretation of our action is as some kind of effective action. In 
classical physics one typically starts from some action principle, obtains the corresponding 
field equations and finally solves them for some initial or boundary conditions. In this 
work we have proceeded in the opposite way: we started from solutions obtained in the 
standard cosmological model and then we have searched for an action that, possessing 
certain properties, gives rise to field equations having the same solutions. 

Classical actions are of course real but effective quantum actions usually have a complex 
structure coming from loops and related to unitarity. The presence of an imaginary part 
in the action, evaluated on some classical configuration, indicates quantum lost of stability 
by particle emission of this configuration [63] . Therefore it is tempting to think that our 
action could have some interpretation in terms of an effective quantum action. However, 
our action determination procedure does not allow to make such a kind of statement. 
Complications in the action could be just an artifact of our construction. In any case it 
was shown that such action exists and it reproduces the present universe evolution without 
DE having = as a vacuum solution. As a drawback of this result, we have shown that 
this f{R) function cannot be considered as a fully consistent gravitational theory since it 
docs not obey the viability conditions revised in Chapter 1 and that therefore it should be 
regarded as an effective model to mimic ACDM cosmological background evolution. 

To conclude the chapter a completely general procedure to reproduce EH gravity with 
an arbitrary perfect fluid by using f{R) theories has been implemented. It has been 
explicitly shown that any perfect fluid when is described by a constant equation of state 
can be mimicked by an appropriate f{R) model. Standard cases for perfect fluids such as 
dust matter, radiation and cosmological constant have been presented in this analysis. 



Chapter 3 

Cosmological perturbations in f{R) 
theories 

3.1 Introduction 



This chapter will be devoted to a study of the evolution of scalar cosmological perturba- 
tions in f{R) theories. To do so, a completely general procedure will be implemented and 
several consequences will be analyzed. The importance for addressing the present problem 
lies in the necessity to discriminate among different DE models, including f{R) modified 
gravities, by using observations. It is well-known that by choosing adequate f{R) func- 
tions, one can mimic any expansion history, and in particular that of the ACDM model. 
Accordingly, the exclusive use of observations from SNIa [64] , baryon acoustic oscillations 
[65] or CMB shift factor [11], based upon different distance measurements which are sen- 
sitive only to the cosmological expansion history, cannot settle the question of the DE 
nature [66]. 

However, there exists a different type of observations which are sensitive, not only to the 
cosmological expansion history, but also to the evolution of matter density perturbations. 
Pioneering work on density perturbations in FLRW cosmological models was presented in 
[67]. The fact that the evolution of perturbations depends on the specific gravity model, 
i.e., it differs in general from that of Einstein's gravity even though the background evo- 
lution is the same, means that this kind of observations will help distinguishing between 
different models capable to explain cosmic acceleration. This is therefore the aim of the 
present chapter: to provide a completely exact method to determine how cosmological 
perturbations grow in f{R) theories and to settle the extracted consequences from this 
result. Such a problem has been exhaustively considered in the literature [68]. In this chap- 
ter we shall show that for f{R) theories the differential equation for the matter density 
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contrast is a fourth order differential equation in the longitudinal (also called conformal 
or Newtonian) gauge but it reduces to a second order differential equation for sub-Hubble 
modes. 

This general equation will be compared with the standard simplification procedure 
(so-called quasi-static approximation) widely used in the previous literature. This appro- 
ximation considers simplifications in the density equation determination from the very 
beginning. We shall show that for general f{R) functions, the quasi-static approximation 
is not justified. However, for those f{R) adequately describing the present phase of acce- 
lerated expansion and satisfying local gravity tests, it does give a correct description for 
the evolution of perturbations. 

Once the general results are presented, some immediate applications may be imple- 
mented. For instance, our analysis may also be used to settle the validity of some proposed 
f{R) models, by comparing the predicted matter spectra with recent observations of LSS 
[69]. 

The present chapter is organized as follows: in Section 3.2 we briefiy review the theory 
of cosmological perturbations for the standard ACDM model, introducing the gauge- 
invariant variables and revising the well-known results for the EH theory in order to es- 
tabhsh a comparison with f{R) gravities. Next, Section 3.3 will be devoted to thoroughly 
study cosmological perturbations in f{R) theories. The perturbed modified Einstein equa- 
tions are obtained through a completely general procedure for those theories in the Sub- 
section 3.3.1. In Subsection 3.3.2 the density perturbations equation is obtained whereas 
the Subsection 3.3.3 compares the obtained exact results with the ones given by using the 
quasi-static approximation. Finally in this section, in Subsection 3.3.4 we shall study the 
growth of perturbations for some particular /(it!) models. Section 3.4 will then show how 
our previous results may be used to constrain or rule some f{R) models out and finally 
some general conclusions for the presented results will be given in Section 3.5. 

This chapter is based upon the results presented in references [70, 71, 72]. 



3.2 Theory of cosmological perturbations 
3.2.1 Generalities 

To study cosmological perturbations, the 4-dimensional full line element may be decoupled 
into background and perturbed parts as follows 

= 9(0) M.dxMx'^ + Sg.^dx^dx-' (3.1) 
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II, v — 0, 1, 2, 3 with g'(Q)^j, representing the homogeneous FLRW background metric and 
5gfj,i, describing a small perturbation. This perturbation may be split [73] in three different 
types: scalar (S), vector (V) and tensor (T) contributions as follows 

Sg^.. = Sg'. + + 5gl.- (3.2) 

This classification obviously refers to the way that fields included in 5g^i, transform 
under three-space coordinate transformations on a constant time hypersurface. Tensor per- 
turbations produce gravitational waves which do not couple to energy density and pressure 
inhomogeneities and propagate freely. Vector perturbations are dumped with cosmological 
expansion and are therefore negligible today. On the contrary, scalar perturbations may 
lead to growing inhomogeneities which will give rise to the large scale structures and the 
CMB anisotropics which are seen today. Thus let us explicitly implement each type of 
perturbations: 

- Scalar perturbations: the most general metric perturbations of this type are given by 
four scalar functions: (f) , , E and B of the space-time coordinates, as follows 

where i, j = 1,2,3 from now on hold for spatial indices and subindex , i means ordinary 
derivative with respect to i*'' -coordinate. 

- Vector perturbations: these perturbations can be represented by two divergenceless three- 
vectors Fj and Si as follows: 

= -Ar) [_\ ^_ -+\_) 

where divergenceless conditions - Einstein's convention applied - mean 

F\i = S\, = (3.5) 

and shift from upper to lower indices -and viceversa- is performed through the spatial 
part of spatially flat background metric tensor, i.e., Sij and its inverse S^^ . 

- Tensor perturbations: tensor perturbations are given by a symmetric three-tensor hij 
satisfying the following conditions 

h\^0 ; h'{j = (3.6) 

i.e., traceless and transversality conditions respectively, meaning that h^j does not contain 
any piece transforming as scalars nor as vectors. Thus, the metric contribution SgJ^^ is 



(3.3) 



(3.4) 
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simply given by 

The number of independent functions introduced to define Sq/^i, without loss of generality 
is ten: four scalar functions for scalar perturbations, two three-vectors for vector pertur- 
bations with one constraint each and one symmetric three-tensor with four conditions for 
tensor perturbations. This number coincides with the number of independent components 
of 5 Qui, as a 4 X 4 symmetric tensor. 



3.2.2 Gauge-invariant variables and gauge choice 

Metric perturbations, as the ones defined above, are gauge-dependent, i.e. an infinitesi- 
mal coordinates transformation could give rise to two apparently different perturbations 
whereas they indeed represent the same physical perturbation. This is the reason why 
Bardeen introduced [74] gauge-invariant quantities that are explicitly invariant under in- 
finitesimal coordinate transformations. The starting point is to consider infinitesimal co- 
ordinate transformations 

^xi^ ^ xi' + ^'^ix). (3.8) 

It can be proven very easily that in the new coordinates {x^} the metric g,j,v{x) can be 
written as: 

~9>.u{x) = 9^u{x) + Ci,g^,{x) + 0{e) (3.9) 

what proves that two metrics g^j^^ and g^^i, differing on a Lie derivative represent the 
same physical perturbation ^. Let us consider a coordinate transformation given by the 
parameters (^°, ^*) , i.e., 

x' = + e = + (3.10) 

where prime holds for derivative with respect to r , and is decomposed as C = + 
5^^ , i.e., it is given by a solenoidal part, ^ , and an irrotational part according to 

Helmholtz's theorem. Therefore dr , dx* and a(r) can be expressed in terms of the new 

^Lie derivative of a twice covariant tensor g^^ with respect to ^ is given by C^^Qfj^v = —g^^x-u — 
9t^\;n + 9fiv;\^^ where 5^^;a = if Levi-Civita connection is considered. 
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coordinates {x^} as: 

dr = df- ^°'df - dx' 
dx' = dx' - C'df - C,j dx^ 

= dx'-(^t + S''e',j)df-{f,j+S^%kj)dx^ 
air) = a(f)-eV(f). (3.11) 

Therefore if identities (3.11) are applied to the expression (3.1), the obtained metric should 
have the aspect of the original line element provided the involved quantities defined in ex- 
pressions (3.3), (3.4) and (3.7) transform as follows: 



- Scalar perturbations: 

^ = V + 

B = B+e-^' 

E = E-^ (3.12) 
where only scalar contributions and ^ are present. 



- Vector perturbations: 

Si = S^ + t (3.13) 
where only vector contribution is present. 



- Tensor perturbations: 

hij = hij (3.14) 

which turn out to be gauge- invariant. 

From previous results, gauge-invariant quantities can be constructed. For scalar per- 
turbations, since four scalar functions were introduced and two scalar gauge parameters 
(^° and ^) are present, two gauge-invariant quantities could be, for instance: 

^ ^ (j)+-\(B-E')a\ ■ ^ ^ ilj + n(B -E') (3.15) 
a 

where by construction $ = $ and ^ — ^ are known as the Bardeen's potentials in [74]. 
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For vector perturbations one gauge-invariant quantity could be 

Si = Si + F-. (3.16) 

With all previous results in mind, one can choose, i.e. one may specify in which coordinate 
system the scalar perturbations are going to be studied. There exist several possibili- 
ties for the gauge choice. Among them we can mention synchronous and longitudinal (or 
conformal-Newtonian) gauges. 

Synchronous gauge: This gauge is defined by the conditions (f) — B — [67]. However, 
it can be shown that the required synchronous coordinates are not completely fixed since 
a residual coordinate freedom remains, what renders the interpretation of calculations in 
this gauge difficult. 

Longitudinal gauge : This gauge is defined by the conditions B = E = and, in this 
gauge, coordinates are totally fixed since E = determines C, uniquely. Using this re- 
sult, B = allows to fix determines C,^ without any uncertainty. We draw the important 
conclusion that in this gauge (p and ip coincide with the gauge invariant variables (3.15) 
$ and respectively which have a simple physical interpretation as the amplitudes of 
the metric perturbations in the usually so-called conformal-Newtonian coordinate system. 



3.2.3 Equations for cosmological perturbations in EH gravity 

In ACDM model within the metric formalism it is possible to obtain a second order 
differential equation for the growth of matter density perturbation. Let us previously 
define the density contrast S as follows: 




Po Po 

where po holds for the unperturbed mean cosmological energy density for a fluid and p 
for the perturbed energy density of the same cosmological fluid. 

In the following, as was mentioned in the beginning of the chapter, the longitudinal 
gauge will be considered to perform our calculations. Thus, the flat FLRW D — A metric 
tensor with scalar perturbations expressed in this gauge and by using conformal time r is 
written as: 

ds=^ = a^ir) [(1 + 2$)dr2 - (1 - 2^){dr^ + r^dnl)] (3.18) 

where $ = $(r, x) and \E' = \E'(r, x) arc the well-known Bardeen's potentials [74]. From 
this metric, the first order perturbed standard Einstein equations are obtained: 



(3.19) 
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Now that the metric with scalar perturbations is known, SC^ is straightforwardly de- 
termined through Einstein's tensor definition. Therefore the first order perturbed 
standard Einstein equations, i.e., for EH gravity read: 



56-° = (2?/$ + 2^')n = STrGa^^T" 

5G) = {-m'^ - 2U^^ - 2W - 2^" - m-^' - V^D)S} + D, ) = SttGo^^T} 

(3.20) 

where D = ^ — , the prime denotes derivative with respect to conformal time r and 
the subindex , i is the usual derivative with respect to the i*'* -spatial coordinate. 

To study the growth rate of cosmological perturbations we shall consider models with 
conventional hydrodynamical matter described by a perfect fluid energy-momentum tensor 
as given in Section 1.2. It should be reminded that = dx^/ds is the mean 4- velocity 
of the fluid. Unperturbed 4-velocity in FLRW conformal coordinates becomes: 

«f„, = -l_(i,o) (3.21) 

and to flrst order in scalar perturbations, it can be shown that 

u^" ^ a{r)-\l-^,5u') (3.22) 

where 5u'^ can be decomposed as follows: 

5u'^v^ + v,' (3.23) 

with If and v,^ being the solenoidal and irrotational components respectively. Note at 
this stage that only contributes to the vector perturbations but not to the scalar ones 
and V is usually referred to as the potential for velocity perturbations. 

Taking into account the previous digression, the perturbed energy-momentum tensor 
components are proven to be: 

ST'', = -ST\ = - (po + Po)diV (3.24) 
with SP pressure fluctuation. Substituting expressions (3.24) in (3.20), they become 



2 • ■ 2 



A7rGa^6p 

4:nGa^{po + Po)v,i 
-4nGa^SPSi 



(3.25) 
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For i j considered, D, — which in Fourier space means kikjO^ — for any i,j 
values and then D is identically null and thus 

$(r, x) = ^(r, x). (3.26) 

Such a result permits to simplify the previous equations (3.25) to become: 

(a$),. = 4:nGa^{po + Po)v,i 
^" + 3W + (n^ + 2n')^ = AnGa^SP. (3.27) 



At this stage, a short digression about the pressure P dependence may be valuable: the 
pressure is, in principle, a qTiantity depending on energy density and entropy per baryon 
ratio. Thus, a pressure fluctuation SP can be expressed in terms of density and entropy 
perturbations as follows 

J' ^{%)/'-'''^''^{^)/' 

where c| can be understood as the squared sound velocity of the fluid perturbations. 
In a single component perfect fluid with constant equation of state there are no entropy 
perturbations. However, if the perturbation description needs to include more than one 
component, entropy perturbations may be present. 

In the following we shall restrict ourselves to adiabatic perturbations, i.e. SS — and 
therefore 

5P = cl6p (3.29) 

and equation of state for the fluid will be considered constant, i.e. P = up with u 
constant. Thus, perturbed and unperturbed content matter are assumed to have the same 
equation of state, i.e. 5P/5p = c| = -Po/po, where cs — for dust matter adiabatic 
perturbations. 

With the previous assumptions, the equations (3.27) can be combined to obtain the 
growth rate 5 evolution in Fourier space ^. For instance, for dust matter the resulting 
differential equation is 



, k^-6pe-18p' k' + 9pi2p-3n')-k\9p-3n') , 

^ ^^k^^w^rm k^-p{3k-+m-) ^ = ° ^^-^^^ 



^In the rest of the present chapter it must be understood that symbols (5 , $ , ^' and v will hold for 
the Fourier corresponding quantities but the subindex k will be omitted in order to simplify the notation. 
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where p = ATiGpoa'^ ~ —T-L' + 'H? according to the background standard Einstein equa- 
tions, as seen for instance from (2.15) setting Pq = and f{R) constant. We point out 
that in order to obtain the equation (3.30) it is not necessary to calculate the potentials 
$ and ^ explicitly, but algebraic manipulations in the equations (3.27) are enough to get 
this result. 

Some limits can be taken in the previous equation: for instance, it is of particular 
interest to consider those k modes whose wavelength is much smaller than the Hubble 
radius. These modes are known as sub-Hubble modes and are identified by the condition 
k ^ Ti or equivalent ly /cr ^ 1 . In this approximation the equation (3.30) reduces to the 
well-known expression: 

S" + W - AnGpoa^S = 0. (3.31) 

In this regime and at early times, the matter energy density dominates over the cosmo- 
logical constant term and it is easy to show that 6 solutions for (3.31) grow as a(r) . At 
late times (near today) the cosmological constant contribution is not negligible and thus 
the equation (3.31) does not admit any more power-law solutions of the type 6 oc a{Ty 
for some 7 . It is necessary in this case to assume an ansatz for 5 : one which works very 
well is that proposed in references [66] and [75] namely 

(3.32) 

where ^M{ci) = and Hm was defined in expression (2.10). Expression (3.32) turns 

out to fit with high precision the numerical solution for 5 with a constant parameter 
7 = 6/11. 



S{a) 



exp 



{nM{ay - l)dlna 



3.3 Cosmological perturbations in f(R) theories 
3.3.1 Perturbed Einstein equations in f{R) theories 

Using the perturbed metric (3.1) and the perturbed energy-momentum tensor (3.24), the 
first order perturbed equations for f{R) theories in the metric formalism, assuming that 
the background equations given in expression (1.10) in Chapter 1 hold, may be written as: 

(1 + /«)5Ge + ((i?(o)),^ + V^V. - 5i:n)fnR6R + [(5^^") V.V, - 5i:{Sg^^)VaV ^Ur 
- "^ro)(5r2:J - Si:g^^^{6r}j] d,fn = -Si:G5Tj^ (3.33) 

where (i?(o))^^ will denote here the Ricci tensor components corresponding to the un- 
perturbed FLRW metric (2.2) in comoving coordinates whose trace provides the scalar 
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curvature already given in equation (2.8). Note that f{R) derivatives with respect to 
i?(o) have been expressed as usually, i.e., /r = df{R)/dR\R^g^ , /rr = d'^f{R)/dR^\R^g^ 
and again □ = V^V" and V is the usual covariant derivative with respect to the un- 
perturbed FLRW metric. Notice also that unlike the ordinary EH case where Einstein's 
equations are second order, the equations (3.33) constitute a set of fourth order differential 
equations. 

For the linearized modified Einstein equations, the components (00) , (ii) , (Oi) = (iO) 
and (ij) , where i,j — 1, 2, 3 , i ^ j , in Fourier space, read respectively: 

(1 + /r)[-A;2($ + ^) - 3?/($' + ^') + {3n' - en^)^ - 3h'^] 

+ 3/r[^(-3$ + ^) - = 2pS (3.34) 

(1 + /r)[$" + + 3n{^' + *') + Sn'^ + {V! + 2%'^)-^] + fR{3U^ - + 3$') 
+/;j(3$ - *) = 2cl~p5 (3.35) 

(1 + fR)[^' + + + *)] + fR{2^ - *) = -2p(l + cl)v (3.36) 

= ^-/^5R (3.37) 

where = R — i?(o) is given by: 
2 



5i? = - 



a2 



3*" + 6{n' + n^)^ + 3n{^' + 3*') - A;^($ - 2*)] (3.38) 



and SP — CgSp has been again assumed. 

By computing the covariant derivative with respect to the perturbed metric V of the 
perturbed energy-momentum tensor , we find the conservation equations: 

V^T;^ = (3.39) 

which do not depend on f{R) explicitly. To first order, the equations (3.39) read 

3*'(1 + 4) -S' + k'^{l + cl)v = (3.40) 

and 

$ + + v' + Hv{l - 34) = (3.41) 

l + Cg 

for the temporal and spatial components respectively. 

In a dust matter dominated universe, (3.40) and (3.41) can be combined to give 

S" + W + - 3*" - 3U^' = 0. (3.42) 
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In this subsection, we are going to obtain the differential equation obeyed by 5 in a dust 
matter dominated universe when the longitudinal gauge and the metric formalism are 
considered to study first order scalar perturbations of f{R) theories. 

To do so, let us consider equations (3.34) and (3.36) for a dust matter dominated 
universe, and combine them to express the potentials $ and ^' in terms of {$', 5, 5'} 
by means of algebraic manipulations. The resulting expressions are the following 

^ - ^(:^{[3(l + /«M*' + ^0 + /«*' + 2p5](l + /«)(7^-/;,^ 

+ [(l + /^)($' + *') + ^(^'-3*')][(l + /fl)(-^'-3^') + 3/;i^]} (3.43) 

and 

* = ^(:^{[-3(l + /«M*' + ^0-3/A*'-2p5][(l + /«)?^ + 2/^ 

x($' + *') + ^{6' - 3*')][(1 + fR)i-k' + 3.n' - en') - mfn]] (3.44) 

where 

V{H, k) = -6(1 + fRfH' + m[fi + 2(1 + fnYn'] + 3(1 + fR)fn{-2n^ + + U'). 

(3.45) 

The second step will be to derive equations (3.43) and (3.44) with respect to r and thus 
and may be rewritten algebraically in terms of {$", 5, S', S"} . These last results 
can be substituted in equations (3.34) and (3.36) to obtain the potentials $ and \E' just 
in terms of {$", S, 8' , 5"} . So at this stage, let us summarize that we have been able 
to express the following quantities 

$ = $($",^";5,5',5") 
m = ^ {<!>", ^"■S,S',S") 

$' = $'($", 5, 5', 5") 

^' = *'($",*"; 5, 5', (5") (3.46) 

but we do not do here explicitly. By the previous expressions we mean that the functions 
on the l.h.s. depend on the functions inside the parenthesis on the r.h.s. in an algebraic 
way. 

The natural reasoning at this point would be to try to obtain the potentials second 
derivatives in terms of {S, 5' , 5"} by an algebraic process. The chosen equations 
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to do so will be (3.42) and the first derivative of (3.37) with respect to r . In (3.42) it 
is necessary to substitute $ and \E'' by the expressions obtained in (3.46) whereas (3.37) 
first derivative may be sketched as follows 



RR 



l + fR 



5R' + 



fRRfn - IrrO- + /r) 



5R. 



(3.47) 



Before deriving, we are going to substitute that appears on (3.37) by lower derivatives 
potentials {$,^,$',^'}, 5 and its derivatives. To do so we consider (3.34) and (3.36) 
first derivatives with respect to r where the quantity v has been previously substituted 
by its expression in (3.40). Following this process we may express as follows 



m" = ^"($, m'- 5, 5', 5") 



(3.48) 



and now substituting the previous result (3.48) in 5R definition given by expression (3.37), 
we can derive equation (3.37) with respect to r . Solving a two algebraic equations system 
with equations (3.42) and (3.47) and introducing (3.46) we are able to express {$", 
in terms of {S, S', S", S'"} . 



= ^"(5, 5', 5", 5'") ; = ^"(5, 5', 5", 5'"). 



(3.49) 



Thus we substitute the results obtained in (3.49) straightforwardly in (3.46) in order to 
express {$, ^, in terms of {6, d', d", d'"} . With the two potentials and their first 

derivatives as algebraic functions of {S,6', 6", 6"'} , we perform the last step: We consider 
$((5, 5, 5", 5"') and derive it with respect to r . The result should be equal to $'(5, 5', S", 5"') 
so we only need to express together these two results obtaining a fourth order differential 
equation for S . Once this fourth order differential equation has been solved we may go 
backwards and by using the results for S we obtain {$", \1'"} from (3.49) as conformal 
time T functions. Analogously from (3.46) the behaviour of the potentials {$, \1'} and 
their first derivatives could be determined. 

The resulting equation for S can be written as follows: 

/34,/5^" + PsjS"' + (a2,EH + MS" + (ai,EH + MS' + {ao,EH + MS = (3.50) 

where the coefficients (i = 0, ...,4) involve terms with and f'^, i.e. terms 

disappearing if the choice fn constant is made. Equivalently, aj^EH {i — 0, 1, 2) contain 
terms coming from EH term and the linear part in i?(o) of f{R) ■ 

It is very useful to define the parameter e = T-L/k since it will allow us to perform a 
perturbative expansion of the previous coefficients a 's and /3 's in the sub-Hubble limit. 
Other dimensionless parameters which will be used are the following: 

nj{i) Ai) 

'^i^yJTi ^ = 1'2'3 ; /^^^ ^ = 1'2- (3-51) 
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where superindex (i) means i^^ derivative with respect to time r . Expressing now the 
a 's and (3 's coefficients as parameter e expansions, we may write 



3 



Oli,EK = 5^ aj,EH « = 0, 1, 2 

fe=l 

A,/ = E/^f/ ^ = 3,4 
fe=i 

8 

A,/ = EA^/ ^ = 0,1,2 (3.52) 
fc=i 

where two consecutive terms in each series differ in a factor. The expressions for the 
coefficients are too long to be written explicitly. Instead, in the following sections we shall 
show different approximated formulae which are proven to be useful in certain limits. As 

a consistency check, wc find that, both in a matter dominated universe and in ACDM , all 
/3 coefficients are absent since fi and /2 defined by expression (3.51) vanish identically. 
For these cases, equation (3.50) becomes equation (3.31) as expected. For instance, in 
the pure matter dominated case, coefficients k, 's are constant and they take the following 
values Ki — —1/2, k,2 — 1/2, ks — —3/4 and K4 — 3/2. 

Another important feature from our results is that, in general, without imposing <C 
1 , the quotient 

+ (3.53) 

Oi2, EH + P2, / 

is not always equal to . In fact only the quotients 



a 



and '-^ (3.54) 



3(1) ""^ 

\f "2,EH 



are identically equal to H . This last result, namely ct^i-En/ Ci2^E}i ' perfect agreement 
with S' coefficient in expression (3.31) when one is studying sub-Hubble modes in ACDM 
theory, i.e., when f3ij i — 0, ...,4 are not present. 



3.3.3 Evolution of sub-Hubble modes and the quasi-static appro- 
ximation 

We are interested in the possible effects on the growth of density perturbations once they 
enter the Hubble radius in the matter dominated era. In this case T-L <^ k and therefore 
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the sub-Hubble limit e <S 1 can be considered. It can be seen that the and Psj 
coefficients are suppressed by with respect to /32j , A,/ and /3oj , i.e., in this limit the 
equation for perturbations reduces to the following second order expression: 



where we have taken only the leading terms in the e expansion for both a and (3 coeffi- 
cients. 

This expression can be compared with the one usually considered in the literature 
by performing the so-called quasi-static approximation, obtained after performing strong 
simplifications in the perturbed equations - (3.34), (3.35), (3.36), (3.37), (3.40) and (3.41) 
- by neglecting time derivatives of $ and ^ potentials, (see [76]). Thus, for instance in 
references [77] and [78], the quasi-static approximation is given by: 



1 J- 1 ^ Jnn 2 



1 + 3 



0? l + 



1 + /. 



0. 



R 



(3.56) 



This approximation has been nonetheless considered as too aggressive in [79] since 
neglecting time derivatives can remove important information about the evolution of per- 
turbations. 

Note also that in equation (3.55) there exists a difference in a power between those 
terms coming from the /-part and those coming from the EH-part. This result differs 
from that in the quasi-static approximation where difference is in a power according 
to expression (3.56). 

In order to compare the evolution for both equations, we have considered a specific 
function 

/te.t(i?) = -4i?°-'' (3.57) 

where H^ units have been used, which gives rise to a matter era followed by a late time 
accelerated phase with the correct deceleration parameter today. In fact this model belongs 
to Class II jiK) models presented in Section 2.4 and it is therefore cosmologically viable. 

Initial conditions in the matter era were given at redshift z — 485 where the EH- 
part was dominant. Results, for k = 0.2/iMpc~^ are presented in Figure 3.1. It can be 
seen that, as expected, both expressions give rise to the same evolutions at early times 
(large redshifts) where they also agree with the standard ACDM evolution. However, at 
late times the quasi-static approximation fails to describe the evolution of perturbations 
correctly. 
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Figure 3.1: 4 with k = 0.2/iMpc ^ for ftest{R) model and ACDM . Both standard quasi-static 
evolution and equation (3.55) have been plotted in the redshift range from 100 to 0. 



Notice that this f{R) model given by expression (3.57) satisfies all the viability con- 
ditions described in Section 1.6 except for the local gravity tests implemented by the 
condition 4 in that section. As is proven below, it is precisely this last condition namely 

< 1 (3.58) 

which in fact ensures the validity of the quasi-static approximation. Therefore we shall 
now restrict ourselves to f{R) models satisfying all the viability conditions, including 

In Appendix A.l we have reproduced all the a 's and the first four /3 's coefficients 
for each 5 term in (3.50). When the sub- Hubble modes are studied and the condition 
<S 1 is imposed, it can shown that the dominant contributions are the first four 
^ coefficients of the / -part plus the first a coefficient of the EH-part for each term in 
equation (3.50). Thus in this case, the full differential equation (3.50) can be simplified as 

04(5'" + cj" + C25" + ci5' + co5 = (3.59) 

where coefficients co,i,...,4 are: 



4 



= lim V/Sl^l A; = 3,4 (3.60) 
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once the condition <^ 1 has been imposed on the corresponding a 's and p 's contribu- 
tions. These coefficients Co,...,4 have been exphcitly written in Appendix A. 2. We see that 
indeed in the sub-Hubble hmit the C4 and C3 coefficients are neghgible and the equation 
can be reduced to a second order expression. Moreover, for our approximated expression 
(3.59) it is true that Ci/c2 = H as can be seen in Appendix A. 2 straightforwardly. 

From those expressions in Appendix A. 2, the second order equation for 5 becomes 



5" + n5'--- 



[ a' 


+ 1(1- 




8 2ki— K2 ) 

9 -2+fi;2 ) 




GfRRk^ 
[ a' 


+ 






h 


8 2ki — K2 ) 

9 -2+K2 ) 


'efRRk^ 
0? 






24 -l+Ki \ 

25 -2+K2 ) 




'&fRRk'^ 

a2 


+ 


i 






24 \ 

25 -2+k;2 ) 



(3.61) 



which can also be written as: 

. ( dfRRk^ 



+ 9\2_ 81 , 
^4) 16 ^ 



9 2ki 
2 



2+K2 



v 



+1) - 



25 
4 



^ "-2+K2 



(1 - Ki)n^6 = 0. 



(3.62) 



Note that the quasi-static expression (3.56) is only recovered in the dust matter era (i.e. 
for % = 2/t) or for a pure ACDM evolution for the background dynamics. Nonetheless, 
in the considered limit <C 1 it can be proven, using the background equations of 
motion, that 



1 -|- «;i — «;2 ~ 



(3.63) 



and therefore 2ki — ^2 ~ — 2 -|- ^2 ~ — 1 -l- ki what allows to simplify the expression (3.62) 
to become (3.56). This is nothing but the fact that for viable models the background 
evolution resembles that of ACDM [31]. 

In other words, although for general f{R) functions the quasi-static approximation 
is not justified, for those viable f{R) functions describing the present phase of accelera- 
ted expansion and satisfying local gravity tests, it does give a correct description for the 
evolution of perturbations. This result has been here stated for the first time shedding 
some light about the controversy which remained about the validity of the quasi-static 
approximation. 



3.3.4 Some proposed models 

In order to illustrate the results obtained in the previous section, we propose two particular 
f{R) theories which allow us to determine - at least numerically - all the quantities involved 
in the calculations and therefore to obtain solutions for the equation (3.61). 
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Figure 3.2: 5k with k = 1.67 /iMpc"^ for f{R) model A evolving according to equation (3.61), 
ACDM and quasi-static approximation given by equation (3.56) in the redshift range from 100 to 0. The 
quasi-static evolution is indistinguishable from that coming from the equation (3.61), but diverges from 
ACDM behaviour as z decreases. 



As was proven before, the background evolution for viable f{R) models resembles 

that of ACDM at low redshifts and that of a matter dominated universe at high redshifts. 
Nevertheless, the f{R) contribution gives the dominant contribution to the gravitational 
action (1.3) for small curvatures and therefore it may explain the cosmological accelera- 
tion. For the sake of concreteness, the models parameters have been fixed by imposing a 
deceleration parameter today qo ~ —0.6 . 

Thus, our first model A will be: 

f{R) = CiRP. (3.64) 

According to the results presented in references [55] and [80] viable models of this type 
belong to Class II introduced in Chapter 2. As was mentioned there, they both include 
matter dominated and late time accelerated eras provided the parameters satisfy Ci < 
and < p < 1 . We have chosen ci = —4.3 and p — 0.01 in units. This choice does 
verify all the viability conditions, including \fji\ <^ 1 today. 



As a second model B we have chosen: 

f(R) = — (3.65) 

with values Ci = 2.5 • 10~^ , ei = 0.3 and C2 = —0.22 also in the same units. 

For each model, we compare our result (3.61) with the standard ACDM evolution and 
the quasi-static approximation (3.56) by plotting S evolution in Figures 3.2 and 3.3 for 
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Figure 3.3: 5k with k = 1.67 /iMpc^^ for f{R) model B evolving according to equation (3.61), 
ACDM and quasi-static evolution given by equation (3.56) in the redshift range from 100 to 0. The 
quasi-static evolution is indistinguishable from that coming from the equation (3.61), but diverges from 
ACDM behaviour as z decreases. 



models A and B respectively. In both cases, the initial conditions are given at redshift 
z — 1000 where 5 is assumed to behave as in a matter dominated universe, i.e. ^^(t) oc 
a(r) with no A; -dependence. We see that for both models, the quasi-static approximation 
gives a correct description for the evolution which clearly deviates from the ACDM case. 

In Figure 3.4 the density contrast evaluated today was plotted as a function of k for 
both models. The growing dependence of 6 with respect to k is verified. This modified 
A: -dependence with respect to the standard ACDM model could give rise to observable 
consequences in the matter power spectrum, as shown in [56, 72], and could be used to 
constrain or even to discard f{R) theories for cosmic acceleration as will be done in the 
next section. 



3.4 A viable f{R) model different from ACDM? 

Some modified f{R) gravity models have recently been proposed (see for instance [58]) 
claiming to be cosmologically viable in spite of having a cosmological behaviour clearly 
distinguishable from ACDM . Contrary to already mentioned opinions which consider that 
self-consistent f{R) gravity models distinct from ACDM are almost ruled out, authors 
in [58] seem to claim that their proposed model would be cosmologically viable. We have 
shown [72] that although that model does satisfy some consistency conditions, precisely 
because of its departure from ACDM behaviour, it does not satisfy local gravity constraints 
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Figure 3.4: Scale dependence of evaluated today = 0) for k/Ho in the range from 10^ to 4-10^. 



and, in acidition, the predicted matter power spectrum conflicts with SDSS data provided 
in reference [69]. The proposed f{R) model reads 



This model does satisfy three of the usual viability conditions for f{R) theories provided 
in [31] and specified in Section 1.6. However, the model (3.66) does not satisfy the fourth 
of those conditions, namely, <C 1 at recent epochs, imposed by local gravity tests 
[33] from solar system. Although it is still not clear what is the actual limit on this 
parameter, certain estimations give < 10~^ today. This condition also ensures that 
the cosmological evolution at late times resembles that of ACDM . However, for the model 
(3.66), Ifnl ~ 0.2 today for a = 2 and ^ 0.25. 

If we arc only interested in considering large scales, local gravity inconsistencies could 
be ignored, but still deviations from ACDM can have drastic cosmological consequences 
on the evolution of density perturbations, as discussed by several authors [56, 70, 79]. 

Thus, the linear evolution of matter density perturbations for sub-Hubble {k ^ Ti) 
modes in ACDM is given by equation (3.31). Notice that in this equation the evolution 
of the Fourier modes does not depend on k . This means that once the density contrast 
starts growing after matter-radiation equality, the mode evolution only changes the overall 
normalization of the matter power-spectrum, but not its shape. 

However, in f{R) theories for sub-Hubble modes, as was thoroughly studied in the 
previous section, the corresponding equation reads as (3.55). Notice from this equation 
the dependence in the 5 term which appears due to the fact that /rr ^ . Moreover, 




(3.66) 



64 



CHAPTER 3. COSMOLOGICAL PERTURBATIONS IN F{R) THEORIES 



Redshift 


EH term (in Ho units) 


/ term (in Ho units) 


RfRR 


100 


61.35 


1.90 • 10-* 


1.45 • 10-*^ 


50 


30.97 


5.34 


1.12 • 10"^ 


20 


12.69 


2.86 • 10*^ 


1.58 • 10-"^ 





-0.45 


1.01 • 10^^ 


0.079 



Table 3.1: Values (in Ho units) for botli the first EH term and the first / term (the one proportional to 
) in the numerator of the 5 coefficient in equation (3.55) for this model (3.66). Different redshifts have 
been considered and the studied scale was k = 0.33 /iMpc"^ . The EH term, which is k independent, 
cannot be ignored at high redshift. In fact, the redshift at which EH and / terms contributions (for 
this value of k ) equal is around z = 45 . The strong suppression with the redshift of the / -part term 
- which is proportional to fc* in the S coefficient in equation (3.55) - comes from the rapid suppression 
of factor as the redshift increases. The dimensionless quantity R fuR has also been included in the 
last column of the table and it is seen to grow as the redshift decreases ( z — )• ) . 

a careful calculation of involved contributions in equation (3.55) shows that in the 5 term^ 
there exist two contributions: one term is proportional to kP' , which is coming from the (3 's 
contribution, i.e. /-part, and another term coming from a 's contribution, i.e. EH-part. 
At high enough redshift the EH-part term dominates, whereas at low redshift the situation 
is reversed and the /-part term becomes dominant. This is the crucial point that explains 
why k -independent terms both in numerator and denominator cannot be straightforwardly 
removed from S coefficient in equation (3.55) but they have to be preserved for a correct 
sub-Hubble modes study. Thus for instance, for k = 0.33 h Mpc"^ , we give explicit values 
in Table 3.1 for the terms in the numerator of the S coefficient in equation (3.55). 

f(R) 

As a consequence the matter power-spectrum P^^ is further processed after equality 
and would differ today from the standard ACDM power spectrum P^cdm xhese two 
quantities would be related by a linear transfer function T{k) given by: 

p/(^) = T{k) Pfc^™ . (3.67) 

This fact changes the shape of the matter power spectrum dramatically, as shown in Figure 
3.5, where normalization to WMAP3 [81] was imposed. In this figure, SDSS data from 
luminous red galaxies [69] and the ACDM power spectrum from the linear perturbation 
theory with WMAP3 cosmological data [81] are also shown. Notice that ACDM gives 
an excellent fit to data with = 11-2, whereas for the f{R) theory = 178.9, i.e. 
13o" out. Even if the overall normalization is drastically reduced by a 20%, which is 
the present uncertainty over this parameter, the discrepancy would still remain at the 7 a 
level. Actually, leaving the power spectrum normalization as a free parameter, the best 
fit would require a 32% normalization reduction and still would be 4.8(7 away as seen in 
Figure 3.5. 



^Let us consider the numerator of this term for simplicity. This does not mean any loss of generality 
since the denominator behaviour is completely analogous. 
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Figure 3.5: Linear matter power-spectra for ACDM and f{R) in [58] with a = 2. Data were taken 
from SDSS [69]. 



The linear transfer function (3.67) for this model has been plotted in Figure 3.6 and 
it has been seen that T{k) follows with a nice fit a power law in the plotted interval 
T{k) = (k/keq)'^-^^ where k^q ^ 10~^ /i Mpc~^ corresponds to the physical scale entering 
the Hubble radius when matter-radiation equality happened. 
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Figure 3.6: Transfer function T(fc) for f{R) model given in equation (3.66) with a = 2 and f^M ~ 
0.25 . The k dependence for T{k) with this parameters choice was seen to be T{k) oc in the plotted 

range of scakis. Both the numerical T(fc) and the power law proportional to k^'^^ were plotted in scales 
k ~ 0.01 — 0.10 /iMpc^^ and no difference was observed. 



3.5 Conclusions 

In this chapter we have studied the evolution of matter density perturbations in f{R) 
theories of gravity. Thus we have presented a completely general procedure to obtain the 
exact fourth order differential equation for the evolution of scalar perturbations in the 
longitudinal gauge. This expression is valid for any general f[R) theory and applicable at 
any scale k . If the EH gravitational action - both with and without cosmological constant 
- is considered in this general expression, well-known standard results are recovered. 

We have also shown that for sub-Hubble modes, the obtained expression reduces to a 
second order differential equation. Hence, we have been able to compare this result with 
that obtained within the quasi-static approximation, widely used in the literature. Our 
research has explicitly shown that for arbitrary /(i?) functions such an approximation is 
not justified. 

However, if we limit ourselves to /(-R) theories for which <^ 1 today, then the 
pcrturbativc calculation for sub-Hubble modes requires to take into account, not only the 
first terms, but also higher-order terms in the e = l-i/k parameter. In that case, the 
resummation of such terms modifies the equation. Thus, this equation can be seen to be 
equivalent to the quasi-static case but only if the universe expands approximately as in a 
matter dominated phase or in a ACDM model. Finally, the fact that for f{R) models 
with <^ 1 the background behaves today precisely as that of ACDM makes the 
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quasi-static approximation correct in those cases. 

We have finaUy apphed our results to prove that no significant departure of f{R) 
theories from ACDM is allowed for those models that intend to be cosmologically viable 
according to recent data. In fact, the strong -dependence appearing in the evolution of 
perturbations has allowed to rule out f{R) gravities which have recently been claimed to 
be cosmologically viable. 
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Chapter 4 



Black holes in f{R) theories 



4.1 Introduction 



We finish the exposition of our research on f{R) modified gravities by considering some 
aspects derived from the study of BHs in these theories. 

Following the motivation already explained in Section 1.1, f{R) models may present 
BH solutions as GR and other alternative gravity theories do. Therefore it is quite natural 
to ask about BHs features in those gravitational theories since, on the one hand, some BHs 
signatures may be peculiar to Einstein's gravity and others may be robust features of all 
generally covariant theories of gravity. On the other hand, the results obtained may lead 
to rule out some models which will be in disagreement with expected physical results. For 
those purposes, research on thermodynamical quantities of BHs is of particular interest. 

These attempts to detect particular signatures from these objects could be experimen- 
tally implemented, as was explained in Section 1.10, at the LHC in the coming years. 
Therefore the generation of mini BHs could provide important information about the co- 
rrect underlying gravity theory. 

Previous literature on f{R) theories [82] proved, by previously performing a conformal 
transformation in the gravitational action, that Schwarzschild solution is the only static 
spherically symmetric solution for an action of the form R + aR^ in D = A. Also by using 
this conformal transformation, uniqueness theorems of spherically symmetric solutions for 
general polynomial actions in arbitrary dimensions were proposed in [83] (see also [84] for 
additional results and [85] for spherical solutions with sources). 

Using the Euchdean action method [86, 87] in order to determine different thermody- 
namical quantities, anti-de Sitter (AdS) BHs in f{R) models have been studied [88]. In 



69 



70 



CHAPTER 4. BLACK HOLES IN F{R) THEORIES 



[89] the entropy of Schwarzschild-de Sitter BHs was calculated for some particular cosmo- 
logically viable models in vacuum and their cosmological stability was discussed. 

BH properties have also been widely studied in other modified gravity theories: for 
instance, [90, 91] studied BHs in Einstein's theory with a Gauss-Bonnet term plus a cos- 
mological constant. Different results were found depending on the dimension D and the 
sign of the constant horizon curvature k. For k = 0,-1, the Gauss-Bonnet term does 
not modify AdS BHs thermodynamics at all (only the horizon position is modified with 
respect to the EH theory) and BHs are not only locally thermodynamically stable, but also 
globally preferred. Nevertheless, for k = +1 and D — 5 (for D > 6 thermodynamics is 
again essentially that for AdS BH) there exist some features not present in the absence 
of Gauss-Bonnet term. Gauss-Bonnet and/or Riemann squared interaction terms were 
studied in [92], where the authors concluded that in this case phase transitions may occur 
with k — —1 . 

Another approach is given by Lovelock gravities, which are free of ghosts and where 
the field equations contain no more than second derivatives of the metric. These theories 
were for instance studied in [93] and the corresponding entropy was evaluated. 

The layout of this chapter is as follows: In Section 4.2 some generalities about BHs in 
f{R) theories such as several aspects of constant curvature solutions for static spherically 
symmetric cases with and without electric charge are presented. Then in Section 4.3 a 
perturbative approach around the EH action, with no previous imposition of constant cur- 
vature, is performed. There we shall find that up to second order in perturbations only BH 
solutions of the Schwarzschild- AdS type are present. Explicit expressions for the effective 
cosmological constant are obtained in terms of the f{R) function. To deal with such 
differential equations we have used the algebraic manipulation package Mathematica [94]. 

Finally, we shall consider in Section 4.4 the BHs thermodynamics in AdS space-time. 
There it will be found that this kind of solutions can only exist provided the theory satisfies 
Rq + f{Ro) < where Rq holds for the constant curvature solution. Interestingly, this 
expression is proved to be related to the condition which guarantees the positivity of 
the effective Newton's constant in this type of theories. In addition, it also ensures that 
the thermodynamical properties in f{R) gravities are qualitatively similar to those of 
standard GR. Then, some consequences for local and global stability for some particular 
f{R) models will be provided in Section 4.5 and figures of thermodynamical regions will 
be shown in Section 4.6. The present chapter is finished by Section 4.7 remarking some 
conclusions of the presented results. 

This chapter is mainly based upon the results that have been presented in references 
[26] and [95]. 
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4.2 Constant curvature black-hole solutions 



The most general static and spherically symmetric D > A dimensional metric can be 
written as (see [96]): 

ds^ = e-^^^''>A{r)dt^ - A-\r)dr^ - r^dQl.^ (4.1) 

or alternatively 

ds^ = A(r)dt' - ii-\r)dr'' - rMQ^_2 (4.2) 

where dO|,_2 is the metric on the S^'"^ sphere. The identification A(r) = e~^*''^M(r) 
and n{r) = A{r) can be straightforwardly estabhshed. For obvious reasons, the $(r) 
function is called the anomalous redshift: notice that a photon emitted at r with proper 
frequency ojq is measured at infinity with frequency oj^o — oJoe~^^'^^ ^/ A{r) . 

Since the metric is static, the scalar curvature R in D dimensions depends only on r 
and it is given, for the metric parametrization (4.1), by: 

R{r) = \{{D-2)[{D-3){l-A{r)) + 2r{A{r)<^'{r)-A'{r))] 

+ [3A'(r)$'(r) - A"{r) - 2A(r)($''(r) - $"(r))] } (4.3) 

where the prime denotes derivative with respect to r . At this stage it is interesting to 
ask about which are the most general static and spherically symmetric metric tensors with 
constant scalar curvature Rq . This metric tensor can be found by solving the equation 
R{r) = Rq . Thus, it is immediate to see that for a $(r) = $o constant, this equation 
becomes 

R{r) = Ro ^ \ [(^2 -5D + 6)(1 - A{r)) + rA'{r){-2D + 4) - r'^A"{r)] (4.4) 



r 



whose general solution is 



A{r) = 1 + a,r'-^ + a,r'-^ - -^^^f" (4-5) 

with oi and 02 arbitrary integration constants. In fact, for the particular case D — A , 
Rq — and $0 = , the metric can be written exclusively in terms of the function: 

= 1 + 7 + ^- (4-6) 

By establishing the identifications ai = —2GM and 02 = , this solution corresponds 
to a Reissner- Nordstrom solution, i.e. a charged massive BH solution with mass M and 
charge Q . Further comments about this result will be made at the end of the section. 
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Now that the general static and spherically symmetric solution for constant curvature 
has been obtained as given by expression (4.5), let us insert the metric (4.1) into the 
general f{R) gravitational action (1.3), and let us also perform variations with respect 
to the functions A{r) and $(r) , in order to find the corresponding modified Einstein 
equations for this parametrization. Thus, we obtain 



and 



(2 -D){1 + /^)$'(r) - r [Wi?'(r)2 + ($'(r)i?'(r) + R"{r))] = (4.7) 



2rA{r)fRRRR'{rf + fnn [2{D - 2)A{r)R'{r) + 2rA{r)R"{r) + A'{r)rR'{r)] 

+ (1 + fn) [2r{A{r)^"(r) - A(r)$'(r)2) + 2{D - 2)A{r)^'{r) - rA"{r) 

+ A'{r){2-D + 3r^'{r))]-r{R + f{R)) = (4.8) 

where as in previous chapters /r = df{R)/dR, /rr = d^f{R)/dR'^ and /rrr = 
d^f{R)/dR^. The above equations look in principle quite difficult to solve. For this 
reason we shall firstly consider the case with constant scalar curvature R — Rq solutions. 
In this simple case the two previous equations reduce to: 



(2-D)(l + fR)^'(r)^0 



(4.9) 



and 



i? + /(i?) + (! + /«) 
+2A{r) ($'(r)2 - $"(r))] = 



A"{r) + ^ — ^ {A'{r) - 2^(r)$'(r)) - 3A'(r)$'(r) 



(4.10) 



As was commented in the Section 1.5, the constant curvature solutions in vacuum of f{R) 
gravities are given by the equation (1.33) which can be rewritten as: 



Rq 



DfiRo 



2(1 + fR{Ro)) - D 



(4.11) 



whenever 2(1 + /7j(i?o)) 7^ D . Thus from (4.9) one obtains $'(r) = and then, for a 
constant scalar curvature Rq , the equation (4.10) becomes 



Ro + f{Ro) + {i + fR{Ro)) 



A"{r) + {D-2) 



A'{r) 



0. 



(4.12) 



By inserting expression (4.11) in the previous equation (4.12), we get 



A"{r) + {D-2) 



A'(r) 



2 

-^Ro 



(4.13) 
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This is a f{R) -independent linear second order inhomogeneous differential equation which 
can be easily integrated to give the general solution: 

A{r) = ci + c^r'-^ - (4-14) 

which depends on two arbitrary constants c\ and C2 . However, this solution has no 
constant curvature in the general case since, as was found above, the constant curvature 
requirement demands Ci = 1 . Then for negative Rq this solution is basically the D 
dimensional generalization studied by Witten [86] of the BH in AdS space-time solution 
considered by Hawking and Page [87]. With the natural choice $o = , that solution can 
be written as: 

= + (4-15) 



where 



with 

2-K— 

being the area of the D — 2 sphere, 

. (4.18) 

is the asymptotic AdS space scale squared and M is the mass parameter usually found 
in the literature. 

Thus we have concluded that the only static and spherically symmetric vacuum solu- 
tions with constant curvature of any f{R) gravity ( it!o < provided) is just the Hawking- 
Page BH in AdS space. However, this kind of sohition is not the more general static and 
spherically symmetric metric with constant curvature as can be seen by comparison with 
the solutions found in expression (4.5). Therefore we have to conclude that there are cons- 
tant curvature BH solutions that cannot be obtained as vacuum solutions of any f{R) 
theory. 

Let us now consider the case of charged BHs in f{R) theories. For the sake of simplicity 
we shall limit ourselves to the D — A case. The action of the theory will be now: 

-5/(R).Ma.well = / d'^A^ + /(^) - ^1^-^""^) (4-19) 

which is a generalization of the Einstein-Maxwell theory. The tensor F^j, is defined as: 

F^, = d^A,-d,A^. (4.20) 
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Considering an electromagnetic potential of the form: = ^y(r), 0^ and the static 

spherically symmetric metric (4.2), we find that the solution with constant curvature Rq 
reads: 

V{r) = 

A(r) = ;x(r) = l-^— + — r . (4.21) 

Notice that unlike the EH case, the contribution of the BH charge to the metric tensor is 
corrected by a (1 + /j?(-Ro)) factor. 



4.3 Perturbative results 

In the previous section we have considered static spherically symmetric solutions with 
constant curvature. In EH theory it is trivial to show that the only static and spherically 
symmetric solution possess constant scalar curvature. However, it is not guaranteed this 
to be the case in f{R) theories too. The problem of finding the general static spherically 
symmetric solution in arbitrary f{R) theories without imposing the constant curvature 
condition is in principle quite complicated. For that reason, we shall present in this section 
a perturbative analysis of the problem, assuming that the modified action given by the 
expression (1.3) is a small perturbation around EH theory. 

The computation we are about to sketch is quite an involved one since, in order 
to calculate the solutions to a given order, it requires to introduce previous order re- 
sults. To deal with such large equations we have used the algebraic manipulation package 
Mathematica [94]. 

Therefore let us consider a f{R) function of the form 

f{R) ^-{D-2)AD + ag{R) (4.22) 

where a <^ 1 is a dimensionless parameter and g{R) is assumed to be analytic when 
expanded for perturbative solutions. Note that nonanalytic functions in a are therefore 
excluded of this analysis. 

By using the metric parametrization given by (4.2) the equations of motion become: 



A(r)(l + fn) {2/i(r) [{D - 2)A'(r) + rX"{r)] + rX'{r)pi'{r)} 

- 2A(r)2 {2/x(r)[(L> - 2)R'{r)fRR + rfuRR R'{rf + rR"{r)fRR] + rR'{r)fj,'{r)fRR} 

- r//(r) A'(r)2(l + /«) + 2rX{ry{R + f{R)) = (4.23) 
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and 

/i(r)|2A(r)i?'(r) [2(D - 2)A(r) + rA'(r)] /rr + r(l + fR){\'irf - 2A(r)A"(r))} 
-A(r)/x'(r) [2(D - 2)A(r) + rA'(r)] (1 + fn) - 2r\{rf{R + f{R)) = (4.24) 

where R = R{r) is given by expression (4.3). Now the A(r) and /x(r) functions appearing 
in the metric (4.2) can be written as follows 

oo 

A(r) = Ao(r) + ^a^A,(r) 

i=l 

oo 

/x(r) = /Xo(r) + ^a>i(r). (4.25) 

i=l 

Notice that Ao(r) and /io('") are the unperturbed solutions for the EH action with cos- 
mological constant given by 

MO = i + ;:d^-(^^ 

Ao(r) = C2/io(r) (4.26) 

which are the standard BH solutions in a D dimensional AdS spacetime (A/) < pro- 
vided). Note that the factor C2 can be chosen by performing a coordinate t reparametriza- 
tion so that both functions, /xo(r) and Ao(r) , could be identified. For the moment, we 
shall keep the background solutions as given in (4.26) and we shall discuss the possibility 
of getting A(r) = //(r) in the perturbative expansion later on. 

By inserting the expressions (4.22) and (4.25) in the equations (4.23) and (4.24) we 
obtain the following 

First order equations: 

[D - 3)//i(r) + r//i(r) + ^— ^ r = (4.27) 



C2 [Ci{D - l)r3-^ - Aor^ + L> - l] g{Ro)r^ + 

+ C2(L>-2)(L>-l)(Ai5r2-L» + 3)//i(r) 

A^r^ 



2A 



D-l 



D 2 

r 



Ai(r) 



+ 



1 + Cir^-^ - 



L>- 1 



[2C2(l-D)r2Ac^K(i?o)+rA;(r)] =0 



(4.28) 
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where gR{Ro) = dg{R)/(i-R\R=Ro and whose solutions are: 

{C,C,-C2Cs)iD-2){D-l) 



Ai(r) = C4D-l){D-2) + 



j.D-3 



- [C,{D - 2)Ad + C2 {g{Ro) - 2A^^^(i?o))] 



C3 , {g{Ro) - 2AngR{Ro)) 2 



Up to second order in a , we obtain the following 
Second order equations: 



(4.29) 



(4.31) 



-C,{D-3y-''- 



_ 2 Apr'' 
D-1 



A2(r) + C2(/^ - 2){D - 1) (-Ai^r' + 1^-3) i^^ir) 



- (^Cir^-^ + r - A'2(r) - 0^0,(0 - 2){D - 1) (-A^r^ + - 3) r''^ 

- C2 [{D - l){C,r'-^ + 1) - Anr'] \2Angn{Ror + 9{Ro) (^-£^^£^4^ - 9r{Ro) 



D-2 



4DAlgR{Ro)gRR{Ro) 



D-2 



- C4 [C,{D - l)r3-^ + 2] [2A,,^^(i?o) - g{Ro)] = 

(4.32) 

where gRR{Ro) = d'^g{R) /dR^lR^R^ and whose solutions are: 



A2(r) = C6 + 
+ 



CeC^ + (C3C4 - C2C,){D -2){D- 1) 



^ + (,(i?o) - 2A.,.(i?o)) {c, + C^gRiRo) - 



(4.33) 



/^2(r) 



C5 , {g{Ro) - 2ADgR{Ro)) (2L>A^y^«(i?o) - (£> - 2)^^(i?o)) 



+ 



(D-2)2(D- 1) 



r". (4.34) 
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Further orders in a^'^'- can be obtained by inserting the previous results in the sub- 
sequent differential equations to get {A3^4^...(r), /i.3_4^...(r)} . Nevertheless, these equations 
become increasingly complicated to be explicitly solved. 

Notice that from the obtained results up to second order in a , the corresponding 
metric has constant scalar curvature for any value of the parameters Ci^2,...,6 • As a matter 
of fact, this metric is nothing but the corresponding one to the standard Schwarzschild- 
AdS geometry, and can be easily rewritten in the usual form by making a trivial time 
reparametrization as follows: 

A(r) = X{r)[-C2iD'' + 3D-2) + C^{D''-3D + 2)a + C6a'' + 0{a^)] 

W) = l^(r). (4.35) 

Therefore, at least up to second order, the only static, spherically symmetric solutions 
which are analytical in a are the standard Schwarzschild- AdS space-times. 



4.3.1 General expression to arbitrary order for constant curva- 
ture 

Let us now assume from the very beginning that the solutions for the equations (4.23) and 
(4.24) belong to Schwarzschild- 74ci5' BH type at any order in the a expansion. Thus we 
can write ( J > provided) : 

A(r) = /x(r) = 1 - + Jr^ (4.36) 



as solutions for the modified Einstein equations (4.23) and (4.24) derived when the mo- 
dification in the gravitational Lagrangian is given by expression (4.22). If we expand the 
quantities Rs and J in terms of parameter a we get: 



oo 

""a 

i=l 

A 



oo 



1=1 



where Rs and Cj are arbitrary constants. The Jj coefficients can be determined from 
expression (1.33), which here becomes: 

R-{D- 2)Ad + ag(R) + 2{D - 1) J(l + ag'{R)) = (4.38) 
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with it! = —D{D — 1) J is the obtained result when calculating R with functions given 
by (4.36). Expanding equation (4.38) in powers of a it is possible to find a recurrence 
equation for the Jj coefficients, namely for the Ji (with I > 0) coefficient, we find: 

(2 - D){D - l)Ji + . ,. / • ^ iJirW ■ ■ ■ {Ji-ir-'9^\Ro) 

i=\j cona.l 

+ 2(D-i)Y,Jk '^' E -j-T^- — ^('/ir(^2)-...(j.-.-i)^--V^+^nit:o) = 

k=0 i=0 cond.2 ^ ^ l K L 

(4.39) 

with Ro = -D{D - l)Jo = DAd and g^^{Ro) = S'^g{R)/dR^^\ r=Rq . In the previous 
recurrence relation, the first sum is done under the condition 1 given by: 

i-i i-i 

= i, e NU {0} and ^mim = l-l (4.40) 

m=l m=l 

and the second one under the condition 2: 

l-k-l l-k-1 

^ = e NU{0} and ^ rnim^l-k-1. (4-41) 

m=l m=l 

For instance we have: 

A{g; D, Ad) 



{D-2){D-1) 



J A{g-, D, Ad)[{D - 2)gR{Ro) - 2DADgRR{R^)] 

{D - 2nD - 1) ^^-^^^ 



where A{g ; D, Ad) = giRo) - 2ADgR{R^) 

Now we can consider the possibility of removing Ad from the gravitational Lagrangian 
(4.22) from the very beginning and still getting an AdS BH solution with an effective 
cosmological constant depending on g{R) and its derivatives evaluated at i?o = . In this 
case the results, order by order in a up to order , are: 

Jo(Ad = 0) 
MAd = 0) 

J2(Ad = 0) 



0, 

gjo) 

{D-2){D - 1)' 

^(0)^^(0) (4 43) 

(D-2)(D-iy ^ > 



As we see, in the context of f{R) gravities, it is therefore possible to have a BH in an AdS 
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asymptotic space-time even if the initial cosmological constant A^i vanishes but ag{0) > 
and a small enough as considered in the explained reasoning. 

To end this section we can summarize by saying that in the context of f{R) gravities 
the only spherically symmetric and static solutions in the general case (without imposing 
constant curvature) in perturbation theory up to second order are the standard BHs in 
AdS space. However, the possibility of having static and spherically symmetric solutions 
with nonconstant curvature cannot be excluded in the case of f{R) functions which are 
not analytical in a . 



4.4 Black-hole thermodynamics 

In order to consider the different thermodynamic quantities for the f{R) BHs in AdS 
space-time, we start from the temperature. In principle, there are two different ways of 
introducing this quantity for the kind of solutions that we are considering here. Firstly 
we can use the definition coming from Euclidean quantum gravity [97]. In this case one 
introduces the Euclidean time t — it and the Euclidean metric ds^ is defined as: 

-dsl = -da^ - r^dnl_2 (4-44) 

where: 

dcr" = e-'*('-U(r)dT2 + A-\r)dr\ (4.45) 

The metric corresponds only to the region r > th where th is the outer horizon position 
with A{rH) — . Expanding dcr^ near th we have: 

^ e-2*(^-U'(r^)pdT2 + (4.46) 



(4.47) 



where p — r — th ■ Now we introduce the new coordinates R and defined as: 
so that: 

dfT^ = R?de'^ + dR^. (4.48) 

According to the Euclidean quantum gravity prescription, r coordinate in expression 
(4.46) is included in the interval defined by and — ^/Te- On the other hand, in 
order to avoid conical singularities, 9 must run between and 27r . Thus it is found that 

T^; = -^e-*('-^M'(rH) . (4.49) 
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Another possible definition of temperature was firstly proposed in [98] stating that 
temperature can be given in terms of the horizon gravity /C as: 



where /C is given by: 



T^^^ (4.50) 



/C = lim ^£^. (4.51) 

'■-^^^ V\9tt9rr\ 



Then it is straightforward to find: 

Tk = Te. (4.52) 

Therefore, both definitions give the same result for this kind of metric tensors. Notice also 
that in any case the temperature depends only on the behaviour of the metric near the 
horizon but it is independent from the gravitational action. By this we mean that different 
actions having the same solutions have also the same temperature. This is not the case 
for other thermodynamic quantities as wc shall sec later. Taking into account the results 
in previous sections for Schwarzschild- AdS BHs, we shall concentrate for simplicity only 
on those solutions, i.e. for a metric as (4.2) where A{r) is given by expression (4.15) and 
$ = is adopted as a natural choice. 

Then, both definitions of temperature lead to: 

'^-f - {D-I)r% + {D-3)P- 

Notice that the temperature is a function of rH only, i.e. it depends only on the BH size. 
In the limit vh going to zero the temperature diverges as T I/th and for th going to 
infinity T grows linearly with tr ■ Consequently T has a minimum at: 



I D - ?> 

rm^hl^—^ (4.54) 



corresponding to a temperature: 



The existence of this minimum was estabhshed in [87] for D — A by Hawking and Page 
long time ago and it is well-known. More recently Witten extended this result to higher 
dimensions [86] . This minimum in the temperature is important in order to set the regions 
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with different tfiermodynamic befiaviours and stability properties. For D = 4 , an exact 
solution can be found for rn '■ 

/ I -\ 2/3 

TH = ^ V3 ■ (4-56) 

181/3 /9«s + ^12 + 81^ J 

Thus, in the Rs <^ I limit, we find th — Rs whereas in the opposite case I <^ Rg , we 
get Th — (fRs)^^^ ■ For the particular case D — 5 , th can also be exactly found to be: 



= o h/l + ^-1 (4.57) 




which goes to i?| for Rs <C / and to IRs for / <C Rs . Notice that for any T > Tq , we 
have two possible BH sizes: one corresponding to the small BH phase with th < Tm and 
the other corresponding to the large BH phase with th > Tm ■ 

In order to compute the remaining thermodynamic quantities, the Euclidean action 



Se — 



IGttG 



l^Jd^x^{R + f{R)) (4.58) 



is considered. Extending to the f{R) theories the computation by Hawking and Page 
[87] and Witten [86] , we evaluate the gravitational Lagrangian in the Schwarzschild- AdS 
scalar curvature solution times the difference between the AdS space-time volume minus 
the BH space-time volume. Thus, we may write: 

AS, = - ^° + /W Ay (4.59) 

where Ro = —D{D — 1)//^ and AV is the volume difference between ^4^5' and BH 
solutions, which is given by: 

/5/iD-2 



AV 



Y) i^'^H-' - rD ■ (4.60) 



2{D- 

Notice that from these expressions, it is straightforward to obtain the free energy F since 
ASe = (3F and therefore 

36n{D-l)GD y H H ) H \ ) 

We see that provided —{Rq + /{Rq)) > , which is the usual case in EH gravity, one has 
F > for Th < I and F < for th > I ■ The temperature corresponding to the horizon 
radius th — I will be denoted Ti and it is given by: 
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Notice that for D > 2 we have Tq < 71 . 

On the other hand, the total thermodynamical energy may now be obtained as: 

dASE {Ro + f{Ro))Mf 

^ ^ ^ ^ — w^) — ^ ^ ^ 

where M is the mass defined in (4.16). This is one of the possible definitions for the BH 
energy for f{R) theories, see for instance [99] for a more general discussion. For the EH 
action with nonvanishing cosmological constant one has f{R) = —{D — 2) A/) and then 
it is immediate to find E — M . However, this is not the case for general f{R) actions. 
Notice that positive energy in AdS space-time requires Ro + f{Ro) < . Now the entropy 
S can be obtained from the well-known relation: 

S^I3E- /3F. (4.64) 

Then one gets: 

^_ {Ro + f{Ro)rAD-2{rH}_ ^^g^^ 



8(D-1)Gd 

where AD_2{rH) is the horizon area given by A£)_2{rH) = t^~^A*d-2- Notice that once 
again positive entropy requires Rq + /(i?o) < . For the EH action with nonvanishing 
cosmological constant one has f{R) — —{D — 2)Kd one has Rq -\- f{Ro) — —2(D — 1)/Z^ 
and then the Hawking-Bekenstein result [100] 

S=%^. (4.66) 
is recovered. Finally, one can compute the heat capacity C which can be written as: 

Then it is easy to find 

^ _ -{Rq + fm){D - 2)/x^-2rg-^/^ {D - l)r% + (D - 3)1' 

8Gd{D - 1) {D- l)rl -{D- 3)P ' ^ '"'^^ 



For the already mentioned case of the EH action with nonvanishing cosmological constant 
one finds: 

(i^ - 2)nn-2r^-' jP - IK + iD- 

AGn {D-I)r%-{D-3)P- ^^'^''^ 

In the Schwarzschild limit Z — > oo , this formula gives: 



AGd 
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which is the well-known negative result for standard BHs of this type. In the general case, 
assuming like in the EH case {Rq + f{Ro)) < , one finds C > for > r^o (the large 
BH region) and C < for th < thq (the small BH region) . For r^j ~ r^o ( T close to 
To ) C is divergent. Notice that in EH gravity, C < necessarily implies F > since 
To < Ti . 

In any case, for f{R) theories with Rq + /(-Rq) < , we have found an scenario similar 
to the one described in full detail by Hawking and Page in [87] long time ago for the EH 
case. 

For T < To , the only possible state of thermal equilibrium in an AdS space is pure 
radiation with negative free energy and there is no stable BH solutions. For T > Tq we 
have two possible BH solutions: the small (and light) BH and the large (heavy) BH. The 
small one has negative heat capacity and positive free energy as the standard Schwarzschild 
BH. Therefore this last configuration is unstable under Hawking radiation decay. For the 
large BH we have two possibilities: if Tq < T < Ti then both the heat capacity and 
the free energy are positive and the BH will decay by tunnelling into radiation, but if 
T > Ti then the heat capacity is still positive but the free energy becomes negative. In 
this case the free energy of the heavy BH will be less than that of pure radiation. Then 
pure radiation will tend to tunnel or to collapse to the BH configuration in equilibrium 
with thermal radiation. 

In arbitrary f{R) theories one could in principle consider the possibility of having 
i?o+/(-Ro) > . However, in this case both the energy and the entropy, given by expressions 
(4.63) and (4.65) respectively, would be negative and therefore in such theories the AdS BH 
solutions would be unphysical. Therefore, Rq + f{Ro) < can be regarded as a necessary 
condition for f{R) theories in order to support the existence of AdS BH solutions. Using 
(1.33), this condition implies l + //j(i?o) > . Let us remind that this condition has a clear 
physical interpretation in f{R) gravities already presented as the condition 2 in Section 
1.5. 



4.5 Particular examples 

In this section we are going to consider several f{R) models in order to calculate the heat 
capacity C and the free energy F since, as was explained in the previous section, these 
are the relevant thermodynamical quantities for local and global stability of BHs. For 
these particular models, Rq can be calculated exactly by using the relation (1.33) with 
R = Ro ■ In the following we will fix the D dimensional Schwarzschild radius in expression 
(4.16) as Rg-^ = 2 for simphcity 

The models that we consider in this section have been previously studied in the lite- 
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rature, but attention there was drawn in studying their cosmological viabihty according 
to conditions provided in Section 1.5. Here we draw our attention on thermodynamics for 
Schwarzschild- AdS BH solutions for these f{R) gravities. 

Both free energy and heat capacity signs are studied for the different values that param- 
eters which appear in these f{R) functions take. Once these signs are known, both local 
or global thermodynamical stability can be determined for these f{R) theories following 
the reasoning explained at the end of the previous section. The considered models are: 

Model I: f{R) = a{-R)f^ 

This model belongs to Class II of f{R) models presented in Section 2.4 if parame- 
ters satisfy a < and < f3 < 1 as seen from Table 2.1. and it could be therefore 
cosmologically viable. 

Substituting in expression (1.33) for arbitrary dimensions we get 



R 



1 _ I) _ (i _ |, 



(4.71) 



Since we are only considering nonvanishing curvature solutions, then we find: 

2-D 



Ro — — 



(2p-D)a 



V(/3-i) 

(4.72) 



Since D is assumed to be larger than 2, the condition (2/3 — D)a < provides well 
defined scalar curvatures Rq . Thus, two separated regions have to be studied: Region 1 
{a < 0, ^ > D/2} and Region 2 {a > 0, ^ < D/2} . For this model we also get 

1 + fniRo) = . (4.73) 

Notice that in Region 1 , 1 + /j?(-Ro) > for D > 2 , since in this case /3 > 1 is 
straightforwardly accomplished. In Region 2 , we find that for D > 2 , the requirement 
Ro+f{Ro) < , i.e. I+Jr^Ro) > , fixes /3 < 1 , since this is the most stringent constraint 
over the parameter /3 in this region. Therefore the physical space of parameters in Region 
2 is restricted to be {a > 0, (3 < 1} . 

In Figures 4.1, 4.2 and 4.3 we plot the physical regions in the parameter space {a, (3) 
corresponding to the different signs of (C, F) . 



Modelll: f{R) = -{-R)" e^p{q/ R) - R 



This model may also belong to Class II of f{R) models presented in Section 2.4 if 
q; = 1 as seen from Table 2.1 and could be therefore cosmologically viable. 
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In this vanishing curvature solution appears provided a > 1 . In addition, we 

also have: 

To get it!o < the condition q{2a — D) < must hold and two separated regions will be 
studied: Region 1 {q > 0, a < D/2} and Region 2 {q <0, a> D/2}. 

In Figures 4.4, 4.5 and 4.6 we plot the regions in the parameter space (a, corres- 
ponding to the different signs of (C, F) . 



Model III: f{R) = R [log{aR)Y - R 

This model may also belong to Class II of f{R) models presented in Section 2.4 if 
q > condition is satisfied and could be therefore cosmologically viable. 

A vanishing curvature solution also appears in this model. The nontrivial one is given 

by 

Since Rq has to be negative, a must be negative as well, accomplishing aRo > . If 
g < is considered then, for D > 2 , the expression (4.75) would imply aRo < 1 and 
then, from f{R) expression for this model, an inconsistency would appear since a negative 
number would be powered to a negative q value. Then q > is the only allowed interval 
for this parameter and therefore there exists a unique accessible region for parameters in 
this model: a < and q > . 

In Figures 4.7 and 4.8 we plot the regions in the parameter space {a, q) corresponding 
to the different signs of {C,F) . 



(4.75) 



Model IV: f{R) = -am, (f )" [l + ^ (f )"]-^ 

As was mentioned in Section 2.5, this model was originally proposed in [33] where it 
was considered to satisfy both cosmological and solar system tests without a cosmological 
constant. For this model, n—1 was considered for simplicity. Hence denoting fR{Ro) = e 
we get 

. -^R^^ (4.6) 

and a relation between mi , D and e can be imposed. Therefore this model would only 
depend on two parameters a and (3 . A vanishing curvature solution also appears in this 
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model and two nontrivial curvature solutions are given by: 




a 



D{mi - 2) + 4 ± ^/f^^/mlD^ - 8D + 16 . 



(4.77) 



2p{D - 2) 



The corresponding 1 + /^(i^o) values for (4.77) are 



4(L>-2)2 



(4.78) 



{VmiD'^ - 8D + 16 ± Dy/ff^Y 



where mi > and mi > {8D — 16)/ D"^ are required for Rq solutions to be real. Since 
1 + fniRo) > is also a condition to be fulfilled, that means that sign(i?Q ) = sign(Q;/3) . 

On the one hand, it can be shown that 1 + //i;(-Rg ) is not positive for any allowed value 
of mi and therefore this curvature solution Rq is excluded of our research. On the other 
hand, 1 + /^(i^o") > only requires mi > for dimension D > A and therefore e < 
is required according to (4.76). Therefore only two accessible regions need to be studied: 
Region 1 {a > 0, /3 < 0} and Region 2 , {a < 0, /3 > 0} . 

In Figures 4.9 and 4.10 we plot the thermodynamical regions in the parameter space 
{a,P) for a chosen e = — 10~^ . Note that 1 + Jr^Rq) does not depend either on a nor 
on (3 and that Rq only depends on the quotient a/ (3 for a fixed mi . 

4.6 Figures for thermodynamical regions 

In the following pages we have plotted accessible thermodynamical regions for previously 
proposed f{R) models I-IV. Thermodynamical regions have been plotted using different 
colors: red: {C < , F > 0} , green: {C > , F > 0} and blue: {C > , F < 0} . 

Parameter spaces have been chosen in order to show possible thermodynamical transi- 
tions between regions. 
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Model I D=4, Region 1 
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(a) Modell, £' = 4,Regionl, a<0, /3>2. (b) Model I, £» = 4 , Region 2, a > , /3 < 1 . 

Figure 4.1: Thermodynamical regions in the plane for Model I in D = 4 . Region 

l(left), Region 2 (right). 
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Figure 4.2: Thermodynamical regions in the plane for Model I in D = 5 . Region 

l(left), Region 2 (right). 
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Model I D= 10, Region 1 
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Figure 4.3: Thermodynamical regions in the {a, /3) plane for Model 1 in D = 10 . Region 
l(left), Region 2 (right). 
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(a) Model II, £> = 4 , Region 1, a <2, q>0. (b) Model II, = 4 , Region 2, a > 2 , q < . 



Figure 4.4: Thermodynamical regions in the {a,q) plane for Model 11 in D = 4. Region 
1 (left), Region 2 (right). 
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Model II D=5, Region 1 
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Figure 4.5: Thermo dynamical regions in the (a, q) plane for Model II in D = 5 . Region 
l(left), Region 2 (right). 
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Figure 4.6: Thermodynamical regions in the {a, q) plane for Model II in D = IQ . Region 
l(left), Region 2 (right). 
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Figure 4.7: Thermodynamical regions in the {a,q) plane for Model III in D = 4 (left) 
and D = 5 (right). 
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Figure 4.8: Thermodynamical regions in the {a, q) plane for Model III in D = 10 . 
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4.7 Conclusions 

In this chapter we have considered static spherically symmetric solutions in f{R) theories 
of gravity in arbitrary dimensions. Firstly, we have discussed the constant curvature case 
(including charged BH in 4 dimensions). Then, the general case, without imposing a priori 
the condition of constant curvature, has also been studied. 

Another important result of this chapter was obtained by performing a perturbative 
analysis around the EH case, assuming regular f{R) functions. We have found explicit 
expressions up to second order for the metric coefficients. These coefficients only gave rise 
to constant curvature (Schwarzschild- ^4^5' ) solutions as in the EH case. 

On the other hand, we have also calculated thermodynamical quantities for the AdS 
BHs and considered the issue of the stability of this kind of solutions. We have found 
that the condition for a f{R) theory of gravity to support this kind of BHs is given by 
Ro + f{Ro) < where Rq is the constant curvature of the AdS space-time. This condition 
also implies that the effective Newton's constant is positive and that the graviton does 
not become a ghost. Consequently thermodynamical and cosmological viabilities of f{R) 
theories turned out to be related as we have shown. 

Finally we have considered several explicit examples of f{R) functions and studied the 
parameter regions in which BHs in such theories are locally stable and globally preferred. 
It was found that the qualitative behaviour is the same as in standard EH gravity but the 
thermodynamical regions are modified depending on the parameters values in each case. 



Chapter 5 

Brane-skyrmions and the CMB cold 
spot 

5.1 Introduction 



Although the properties of BH solutions in BW models have been exhaustively studied in 
the literature, this is not the case of a different type of (topological) spherically symmetric 
solutions, the so-called brane-skyrmions. These configurations may appear in a natural 
way in a broad class of BW theories. 

On the one hand, this type of textures can be understood as holes in the brane which 
make possible to pass through them along the extra-dimensional space. For such objects, 
cosmological involved scales and consequences have not been already studied in the litera- 
ture. Therefore, the leitmotiv of this chapter will be precisely to study some astrophysical 
and cosmological effects that those configurations may have in the CMB data. 

On the other hand, some striking features have recently been discovered [101] in the 
CMB data such as the presence of an anomalous cold spot (CS) in the WMAP temperature 
maps. This CS has been interpreted in different manners by making use of different physical 
mechanisms. Among them, one of the most exciting from the physical point of view could 
justify the appearance of cold spots as the result of a collapsing texture, perhaps coming 
from some early universe Grand Unified Theory (GUT) phase transition. 

In this chapter we shall propose an alternative explanation to those in the existing lite- 
rature: it will be shown that the brane-skyrmions provide a natural scenario to reproduce 
the CS features and that typical involved scales, needed for the proposed brane-skyrmions 
to describe correctly the observed CS, can be as low as the electroweak scale. 
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The present chapter is thus organized as follows: brane-skyrmions in brane- worlds are 
presented in Section 5.2 and then we shall determine the possible distortion in the fractional 
profile of temperatures generated by the presence of these objects. Next, in Section 5.3 we 
shall present the most recent results that claim the existence of a CS in the CMB data. 
Different possible explanations that have been previously proposed in the literature will 
be summarized in Section 5.4. 

The subsequent physical interpretation of the calculations performed in Section 5.2 
will be presented in Section 5.5. In this section, we will also show that those results are 
in agreement with the CMB data and with other theoretical proposals. To conclude the 
chapter, possible future detection of these brane-skyrmions and some conclusions will be 
studied in Section 5.6. 

The results of this chapter were originally published in [102]. 



5.2 Spherically symmetric brane-skyrmions 

In this section we are going to generalize the results presented in Section 1.8 for static 
brane-skyrmions. Here we shall allow time dependence for these objects but spherical 
symmetry will be preserved. Starting from the physical branon fields equations provided 
in equation (1.65), one may define the brane-skyrmion spherical coordinates with winding 
number nw by the nontrivial mapping tt" : — > as follows: 

(l>K^cf>, Ok ^9, XK^F{t,r) (5.1) 

with boundary conditions satisfying the requirement 

F{t, oo) - F{t,0) ^nw7^- (5.2) 

This map is usually referred to as the hedgehog ansatz. With these coordinates (5.1) 
introduced in the metric (1.58), the Nambu-Goto action given in expression (1.53) may be 
rewritten in terms of the previous coordinates as follows: 



-f J d'x Sine (2j^sm (F) cos (F)) 1 - ^ (f^ - F") 



1/2 

(5.3) 



Varying this action with respect to the function F{t, r) the equation of motion for the 
skyrmion profile is obtained and it becomes: 

2 F-F" + ^ (fF'^ + F"F^ - 2F'FF') 
sin(2F) - 2rF' + (r^ + ^^sin^F ) r ^ = (5.4) 
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where dot and prime denote - throughout this section - derivatives with respect to t and 
r respectively. 

In this chapter we are interested in the potential cosmological effects due to the presence 
of a brane-skyrmion within our Hubble radius. For that purpose, gravitational field pertur- 
bations will be computed at large distances compared to the size of the extra dimensions, 
i.e., we are interested in the region 

» Rl = ^- (5.5) 

Notice that in this region, gravity behaves essentially as in four dimensional space-time 
and standard GR can be used in the calculations. Notice also that in order to simplify 
these calculations, the effects due to the universe expansion will be ignored or at least 
assumed negligible. This assumption is fully justified provided r <S Hq^ . In such a case 
the unperturbed (ignoring the presence of the defect) background metric can be taken as 
Minkowski, i.e., g^^i, = rj^y in expression (1.58). From assumptions in (5.5), the equation 
of motion (5.4) reduces to: 

(Fo - + sin(2Fo) - 2rF^ = (5.6) 

which is equivalent to expression (3) in reference [103]. Notice that this is an expected 
result since, as shown in [104, 105], at large distances, i.e. except in the microscopic 
unwinding regions, the dynamical evolution of the fields is completely independent from the 

symmetry breaking mechanism, it simply depends on the geometry of the coset manifold 
K . On small scales [104, 105] it is possible that higher- derivative terms could affect the 
dynamics and even stabilize the textures, this is also the case of brane-skyrmions [49], 
although generically they could unwind by means of quantum-mechanical effects. 

The equation (5.6) admits an exact solution with winding number nw equals to unity 
of the following form: 

Fo(i,r) = 2arctan (^-^) (5.7) 

with i < since as explained below, we will be interested in photons passing the texture 
before it collapses. 

Our approximated equation (5.6) is consistent with the complete equation (5.4) for the 
above solution (5.7) since the second term in the numerator in equation (5.4) vanishes for 
F = Fq and 

^0 - = ^^^2 ; = 4 (5.8) 

so that in the considered regime, the neglected terms in equation (5.4) are indeed irrelevant 
for all r and t values. 
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Now that Fo{t, r) has been determined, we may calculate the energy-momentum tensor 
components also in this region from the Nambu-Goto action (1.53) as: 

J.,. = 2|^|-V2^ (5.9) 

where note that g is the involved metric to determine this tensor. In spherical coordinates 
they become 

Tqq — — ^7T^ ; Tri 



Av^rt 2t;2r2(t2 _ ^2) 



2 



(t2+r2)2 ' (t2+r2)2 

T^^ = sm^OTee (5.10) 
and note that V^T^ identically vanishes. 

We shall now determine the background metric g^y in the r ^ Rb region and in 
the presence of the brane-skyrmion as a small perturbation on the Minkowski metric, i.e. 

Thus, for the scalar perturbation of the Minkowski space-time in the longitudinal gauge, 
the line element will adopt the form, straightforwardly obtained from equation (3.18) in 
Chapter 3 if a(T) = 1 , that follows: 

ds^ = (1 + 2^)df - (1 - 2*)(dr2 + rW^) (5.11) 

and thus the perturbed Einstein tensor components in cartesian coordinates are the follow- 
ing (see [106]): 

5G° = 2di^ (5.12) 

with i,j — 1,2,3 and = X]i=i ^i^* ■ Using Einstein's equations 5G^ — —SttGT^ we 
determine that the potentials $ and ^ are 

^ = $ = AnGvHogi ] . (5.13) 

The physical metric on which photons propagate is not the g^i, we have just calculated, 
but the induced metric (1.58). However, using the solution in (5.7), we find that the 
contribution from branons fields is 0{R'^/r'^) , i.e. neghgible in the region we are interested 
in, so that g^^ ~ . 
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Figure 5.1: Schematic representation of photon trajectory passing near the brane-skyrmion 
as originally presented in [103] . it! is the impact parameter (the radius in the plotted circle) 
and the photon trajectory through the bottom horizontal line is chosen along the Z axis. 



Photons propagating on the perturbed metric will suffer red (blue)-shift due to the 
Sachs- Wolfe (SW) effect [106]. The full expression for the temperature fluctuation is given 
by: 



SW 



= -[^n + ^ ' + dr = -[mi + ' 26dT (5.14) 



where we have considered local and integrated SW effects and neglected the Doppler contri- 
bution. Tj and Tf are the initial and final times respectively required to study temperature 
fluctuation in that interval. Substituting expression (5.13) in the previous one and using 
that 

= z'^ + R'^ ; z = t-to (5.15) 

as may be seen at the schematic representation in Figure 5.1, the fractional distortion we 
obtain is: 



( ^] = 8nGv^ ( , ^° arctan ( f ° + \ - log|io + ;^|^ (5.16) 

where R is the impact parameter and to is the time at which the photon passes the 
texture position at z = . In the limit where Zi ^ —oo and Zf ^ oo the result is 




(5.17) 



with e = Sn'^Gv'^ . As shall be seen in the next section, this temperature profile is the one 
which fits the apparently observed anomaly in the CMB data. 
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VJCI.I, scale = 300' 





Figure 5.2: The combined and foreground cleaned Q-V-W WMAP map after convolution 
with the SMHW at scale Rg as shown in [115]. The CS position is marked and seen at 
the bottom left part of the plot. 

5.3 Cold spot in WMAP data 

One of the most important pieces of information about the history and nature of our uni- 
verse comes from the Cosmic Microwave Background (CMB). Measurements of the CMB 
temperature anisotropics obtained by WMAP [11, 107] have been thoroughly studied in 
recent years. Such anisotropics have been found to be Gaussian as expected in many stan- 
dard cosmological scenarios corresponding to density fluctuations of one part in a hundred 
thousand in the early universe. However, by means of a wavelet analysis, an anomalous CS, 
apparently inconsistent with homogeneous Gaussian fluctuations, was found in [101, 108] 
in the southern hemisphere centered at the position b = —57 ° , I = 209 ° in galactic 
coordinates. ^ The characteristic scale in the sky of the CS is about 5 ° . 

The CS position in usual WMAP plots is shown in Figure 5.2. The existence of this 
CS has been claimed to be confirmed more recently in reference [109]. However, in recent 
references [110] it has been argued that there is no compelling evidence for deviations from 
the ACDM model in the WMAP data. In particular, it is claimed that the evidence that 
the CS is statistically anomalous is not robust. 

^where b is the galactic latitude measured from the plane of the galaxy to the object using the Sun as 
vertex. The galactic longitude is referred to as I and it is measured in the plane of the galaxy using an 
axis pointing from the Sun to the galactic center. 
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Concerning the origin of the CS, different explanations have been proposed in the 
hterature: in references [111, 112, 113] the CS was tried to be explained by voids in the 
matter distribution whereas the reference [114] dealt with possible string theory bubble 
origin. Other authors proposed in [115] the Sunyaev-Zeldovich (SZ) effect as possible 
explanation showing that the flat frequency dependence of the CS is incompatible with 
being caused by a SZ signal alone. However, a combination of CMB plus SZ effect may 
explain the spot and could have a sufficiently flat frequency spectrum. 

5.4 Cold spot as a cosmic texture 

A very interesting possibility about the CS origin has recently been proposed in [116]. 
According to it, some theories of high energy physics predict the formation of various 
types of topological defects, including cosmic textures [104] which would generate hot 
and cold spots in the CMB [103]. These textures would be the remnants of a symmetry 
breaking phase transition that took place in the early universe. 

In order to produce textures, the cosmic phase transition must be related to a global 

symmetry breaking pattern from one group G to a subgroup H so that the coset space 
K = G/H has a nontrivial third homotopy group. A typical example is = SU{N) , 
which has associated n3{K) = Z for > 2 . Notice that, as usual, in order to get 
a texture formed in the transition, the symmetry breaking must correspond to a global 
symmetry since if it were local it could be gauged away. 

Textures can be understood as locahzed wrapped field configurations which collapse 
and unwind on progressively larger scales. These textures can produce a concentration 
of energy which gives rise to a time dependent gravitational potential. CMB photons 
traversing the texture region will suffer a red (blue) shift producing a cold (hot) spot 
in the CMB maps. In reference [116] the authors consider a SU{2) non- linear sigma 
model to build up a model of texture that could explain the observed CS. They simulate 
the unwinding texture by using a spherically symmetric scaling solution and they find a 
fractional temperature distortion given by: 



where 6 is the angle from the center, e is a measure of the amplitude and 6c is the scale 
parameter that depends on the time at which the texture unwinds. The best fit of the 
CS is found for e = 7.7 x 10~^ and 6c = 5.1°. Furthermore, the parameter e is given 
by e = Stt^G^o j where $o is the fundamental symmetry breaking scale which is then set 
to be $0 — 8-7 X 10^^ GeV. This scale is nicely close to the GUT scale thus making the 
results given in [116] extremely interesting. 




(5.18) 
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Nonetheless, it is important to stress that textures require having a global symmetry 
breaking but what one finds typically in GUTs is a local symmetry breaking producing 
the Higgs mechanism and then destroying the topological meaning of texture or any other 
possible defect appearing in the cosmic transition. 



5.5 Physical interpretation of the results and involved 

scales 



As has been shown in Sections 5.2 and 5.4, the results for the fractional temperature 
distortion provoked by the CS given by expression (5.18) coincides with the fractional 
distortion caused by a brane-skyrmion given by expression (5.17) provided that: 

*. ^ „ ; 2 (i) ^ 5 (5.19) 

where 9c was defined in reference [116]. 

Let us now estimate the required scales in order for this kind of branc-skyrmion to be 
able to explain the observed CS. For the minimal model supporting brane-skyrmions with 
three extra dimensions 6 = 3, the relation (1.50) approximately becomes 

M| ~ RlMl (5.20) 

since gravitation is embedded in D = 7 dimensions and therefore R% is the characteristic 
volume for the compactified extra dimensions. Therefore 



and consequently 



If the parameter e needs to be fixed around 7.7 • 10~^ , then the required v would be 

v~ 1.2-lO^^GeV (5.23) 

which in fact can be achieved with ~ / ~ TeV. Notice that for this parameter range, 
the radius of the extra dimension is around 

i?B~10-^m, (5.24) 



5.5. PHYSICAL INTERPRETATION OF THE RESULTS AND INVOLVED SCALESlOl 



i.e., previous approximations in expression (5.5) are totally justified. 

In addition to the previous estimations, we have also checked that the possible effects 
coming from a nonvanishing branon mass are suppressed by 0{M^/v^) . Therefore for 
mass values also around M ~ TeV, which are typical of branon as dark matter candidate, 
see [25] , such effects are neghgible. 

In other words, brane-skyrmions provide an accurate description for the CS without 

any need of introducing very high energy (GUTs) scales since the correct temperature 
ffuctuation amplitude can be obtained with natural values around the electroweak scale. 



5.5.1 Brane-skyrmions abundance 

The abundance of brane-skyrmions in this model is given by its low-energy (large distance) 
dynamics as described by a non- linear sigma model as that in [116]. This is nothing but the 
well-known fact that, except in the microscopic unwinding region, the field evolution only 
depends on the geometry of the coset space K , but not on the details of the symmetry 
breaking mechanism. As was shown in [105] in a simple model with a potential term, the 
final abundance of defects and other properties of the pattern of density perturbations are 
expected to be not very sensitive to the short distance physics, once the texture unwinds, 
making this kind of theories highly predictive. For that reason, we expect that provided 
the same kind of initial conditions are imposed in both models, the predicted abundance 
of hot and cold spots agree with that obtained from simulations in [116]. Such simulations 
show that the number of unwinding textures per comoving volume and conformal time r 
can be estimated as 

dr t4 ^ ^ 

with z/ ~ 2 . This allows to estimate the number of hot and cold spots in a given angular 
radius interval. As commented above, this is a quite robust result at late times with little 
effects for short distance dynamics. In the case of brane-skyrmions, the short distance 
effects will be embodied in the higher-derivative terms appearing in the expansion of the 
Nambu-Goto action or even in possible induced curvature terms generated by quantum 
effects [49]. Provided that such terms do not stabilize the brane-skyrmions, we expect that 
such an abundance could be directly applied in our case. 
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5.6 Future prospects and conclusions 

In this section we dedicate some lines to mention how both present and future experiments 

may confirm the vahdity of the BW models. Then we shall sketch the most important 
conclusions and consequences of the brane-skyrmion model presented in this chapter as a 
viable explanation for the observed CS. 

On the one hand, the fact that the fundamental scales of the theory are of the order of 
the TeV opens the possibility to test this explanation with collider experiments through 
the production of real or virtual branons and KK-gravitons. The expected signatures of 
the model [117] from the production of KK-gravitons come fundamentally from the single 
photon channel studied by LEP, which restrict Md > 1-2 TeV at 95% of confidence level. 

On the other hand, the LHC accelerator will be able to test the model up to — 3.7 
TeV, analysing single photon and monojet production to establish the number of extra 
dimensions as was explained in Section 1.10. 

Finally, there also exists the possibility to find signatures of the model at low energies 
associated with branon phenomenology. This case is more interesting from the cosmological 
point of view, since branons can constitute the nonbaryonic dark matter abundance as 
typical WlMPs [25, 118]. Present constraints coming from the single photon analysis 
reahzed by L3 (LEP) imply / > 122 GeV (at 95 % C.L.) [119, 120] and the LHC will be 
able to check this model up to / = 1080 GeV through monojet production [118]. The 
idea to test the physics associated with the CS with the next generation of colliders at the 
TeV scale is a very intriguing and distinctive property of this texture. 



Chapter 6 

Conclusions and prospects 



In this thesis we have studied some cosmological and astrophysical consequences in two di- 
fferent types of modified gravity theories: f{R) theories and brane- world extra dimensions 
theories. 

First, we have shown some relevant results for f{R) modified gravities. Thus, we 
have studied how some features of general relativity can be mimicked by f{R) models if 
certain conditions are imposed on those models. We have also proved that there exists a 
class of f{R) models which can indeed reproduce the present cosmological evolution -from 
matter-radiation equality till today - as described by general relativity within the standard 
cosmological concordance model. These f{R) models allow to remove the cosmological 
constant term and therefore the observed acceleration is understood to have a purely 
geometrical interpretation with no dark energy origin. Furthermore, initial conditions may 
be imposed on such a class of theories in order to recover a null scalar curvature solution 
for a vacuum scenario. Unfortunately such functions arc not cosmologically viable and 
should be considered as effective models to reproduce ACDM evolution. 

Also, it has been proved that any perfect fluid parameterized by a constant equation 
of state can be reproduced by the presence of f{R) terms. Thus, these functions could 
provide a mechanism to reproduce the cosmological behaviour of dark energy type fluids. 
It therefore seems worthwhile to try to search for a unique phenomenological f{R) theory 
able to explain the main features in the cosmological evolution, from inflation till late time 
observed acceleration. 

Next, we have considered the modification introduced by the new f{R) terms in the 
evolution of cosmological perturbations. Here, an analysis concerning first order cosmolog- 
ical scalar perturbations has been performed. Wc have presented a general method in order 
to obtain the evolution equation of the matter density contrast for arbitrary f{R) theories. 
That procedure was proven to be valid regardless of the chosen f{R) model or the size of 
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the involved scales. An important by-product of the calculations is the conclusion that the 
usual approximations made in the literature for sub-Hubble modes are not always valid. 
In fact, the required hypotheses to get the so-called quasi-static approximation have been 
established here. For those sub-Hubble modes, it was found that only for f{R) models 
satisfying the local gravity constraints, the evolution of perturbations is indistinguishable 
from that obtained using the quasi-static approximation but perfectly distinguishable from 
the evolution obtained by the ACDM model. 

Now that those theoretical calculations were made, the correct evolution for matter 
perturbations in sub-Hubble modes for f{R) models is obtained. Therefore some robust- 
ness tests may be used to compare experimental data with expected theoretical results. 
Thus both the validity and certain constraints on those f{R) functions may be estab- 
lished. This procedure shows how our theoretical results may be used in order to confront 
the predictions made by a given f{R) theory with large scale structure observations. 

As a final subject within our f{R) models research, we dealt with some features of black 
holes in arbitrary dimensions. First, the scalar constant curvature case for static and spher- 
ically symmetric solutions in vacuum was studied. It was determined that the only possible 
solution for the modified Einstein equations in this case was the Schwarzschild-(anti-)de 
Sitter solution. As a complementary result to the previous research, we studied the static 
and spherically symmetric case without imposing the constant curvature condition. To 
do so, a perturbative approach was performed around the standard Schwarzschild-anti-de 
Sitter solution of general relativity. It was found that, up to second order in perturbations, 
the only sohition was the generalized Schwarzschild-anti-dc Sitter solution with modified 
coefficients in terms of f{R) and its derivatives evaluated at the background curvature. 

This part of the thesis was finished by studying the thermodynamic properties of 
Schwarzschild-anti-de Sitter black holes in f{R) theories. Thus accessible thcrmodynami- 
cal regions for several f{R) models were studied in detail. In that realm, a very interesting 
property, which encourages further investigation, was found: it turned out that thermo- 
dynamical viability of f{R) theories in constant curvature solutions was embodied by the 
condition 1 -|- df{R)/dR > 0. This condition was also required to ensure gravitational 
viability for f{R) models. Consequently two different aspects of the viability of f{R) 
theories have been brought together. 

The last part of the thesis was devoted to a study of some consequences of the existence 
of topologically nontrivial brane configurations, the so-called brane-skyrmions, in brane- 
world theories. In particular we have considered the effects of brane-skyrmions on the 
temperature fluctuations of the CMB and their potential connection with the so-called 
CMB cold spot. We have shown that brane- world theories can naturally accommodate 
this striking cosmological feature by invoking these nontrivial topological configurations. 
A different type of textures related to Grand Unified Theories has also been considered 
in the literature, but the model presented here shows that this type of defect could be 
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explained from the electroweak scale physics. 

To conclude, the results included in this thesis have shown that modified gravity theo- 
ries remain compelling candidates to describe the properties of the gravitational interaction 
on very large scales. Both the validity and viability of these theories have still to be sub- 
jected to many theoretical and experimental tests. To do so, the possibility of reproducing 
standard cosmological results, first order perturbations, black holes and features on the 
CMB temperature maps were considered as interesting aspects which could shed some light 
on the fundamental properties of the gravitational interaction. 
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Appendix A 



Coefficients in f(R) cosmological perturbations 

In this appendix we include the coefficients for the series expansions given in equation 
(3.50) in Chapter 3. All the a 's terms (coming from EH-part) and the ffist four f3 's 
terms (coming from /-part) for each S derivative, i.e. 5'^'", ...,5, , in equation (3.50) have 
been included once the condition \fji\ <^ 1 has been imposed to simplify their expressions. 
For sub-Hubble scales the excluded /3 's terms are negligible with respect to the written 
ones. 



A.l Appendix I: a's and /S's coefficients 

Coefficients for 5™ term: 

^« SfUl + fRffte' 

/3g ~ 216fRhe%-2 + H2f. (A.l) 

Coefficients for 5"' term: 

fii'j 8fUl + fRfftne'[3 + fRi3 + h)] 

f3i'j ~ 6f^f^'He\8M-2 + K,)+4Ml2K, + 9K,-2{9 + Ks)]} 

/3g ~ -72/2/i?^e6(-2 + «2)[-4/2(-2 + «2) + /i(19-23«i-10«2 + 4«3)] 

45 ~ _216/^7{e«(-2 + «2)'[-2/2(-2 + «:2) + /i(7-ll«:i-4«2 + 2«3)]. (A.2) 

Coefficients for 6" term: 

«Sh = 432(1 + /«)^°?{V(-l + «i)(-2 + /.2)=^ 
«Sh = 1296(1 + fRy'nh''{-l + K,f {-2 + K2f 
«Sh = 3888(1 + /ij)i°?/V2(-l + «i)'(-2 + «2)' 

/3g ^ sfUi + fRfftn' 

/3g ^ 88/l/f?/^6^(-2 + ;.2) 

~ 24fyi'H\\-2 + n2){-28 + 2Ki + 13K2) 
/3g ~ 72/^/i?^V(-2 + «2)'(-14 + 4«:i + 5«;2). (A.3) 
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Coefficients for S' term: 

"Sh = 432(1 + /^yoH-V(-l + Ki)(-2 + «:2)' 

myi'H'e'i-2 + K,) 
/^g ~ 24/2/27{V(-2 + fi;2)(-28 + 2^1 + 13^2) 

~ 72/^/i?^V(-2 + «:2)'(-14 + 4«:i + 5«:2). (A.4) 

Coefficients for 5 term: 

«o!eh = 432(1 + /«)io?{V(-l + «i)(2«i-«2)(-2 + «2)=^ 

= 1296(1 + /r)^°?{V°(-1 + «i)2(-1 + 4«i-«2)(-2 + «2)' 

43,H = 3888(1 + /«)1°?^V2(-1 + «:i)2(2k2-k2)(-2 + «:2)' 

^ -y/l(l + /i^)Vi^l2 + /^(2 + 2A-/2-2«:i)-2/.i] 

45 :^ 112/|/37^V(-l + «i)(-2 + «:2) 

45 ~ 48/^/i27^V(-l + «i)(-2 + «2)(-16 + 2«i + 7«2) 

45 ~ 144/H/i-HV(-l + «:i)(-2 + «;2f(-6 + 4ACi + «;2). (A.5) 



A.2. APPENDIX II ; C'S COEFFICIENTS 
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A. 2 Appendix II: c's coefficients 

The coefficients appearing in equation (3.59) are explicitly provided here once the condition 
Iful <^ 1 has been imposed. 

X [2/fl/2A;' + 3n^{-7 + lUi + 4^2 - 2^3)] + 2fRf{K^k''{-Q + 6ki + 3«2 - ^s)} 

C2 = [-/iiM'-37{2(-2 + «;2)]' 

+ 5fRfiH^k^{-2 + K2) + m\-l + «i)(-2 + «;2)] 

ci = [-/iiM'-37{2(-2 + «;2)]' 

x[flfik^n + 5fRhn^e{-2 + K2) + m\-i + «i)(-2 + «2)] 

Co = ^?{2(-l + «i)[-/R/iA;^-37{2(-2 + «2)]'[2/l/iV + 9/fl/i?{V(-2 + «2) 

+ 9?^'(2/ti-fi:2)(-2 + fi;2)]. (A.6) 
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